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Chapter 1 

Introduction 



Many approaches to three-dimensional constrained macromolecular chains at thermal 
equilibrium, at about room temperatures, are based upon constrained Classical 
Hamiltonian Dynamics (cCHDa). Quantum- mechanical approaches (QMa) have also 
been treated by different researchers for decades. QMa address a fundamental issue 
(constraints versus the uncertainty principle) and are versatile: they also yield classical 
descriptions (which may not coincide with those from cCHDa, although they may 
agree for certain relevant quantities). Open issues include whether QMa have enough 
practical consequences which differ from and/or improve those from cCHDa. We shall 
treat cCHDa briefly and deal with QMa, by outlining old approaches and focusing 
on recent ones. In QMa, we start with Hamiltonians for N(~^> 1) non-relativistic 
quantum particles, interacting among themselves through potentials which include 
strong vibrational ones (constraining bond lengths and bond angles) and other weaker 
interactions. We get (by means of variational calculations) effective three-dimensional 
constrained quantum partition functions at equilibrium (Zq) and Hamiltonians (Hq) 
for single-stranded (ss) macromolecules (freely-jointed, freely-rotating, open or closed) 
and for double-stranded (ds) open macromolecules. Due to crucial cancellations, we 
can neatly separate the constrained degrees of freedom (by getting the large constant 
vibrational zero-point energies associated to them) from the slow unconstrained angular 
variables (accounted for by Zq and Hq). In the classical limit, we obtain classical 
partition functions Zq from Zq. The Zq^s are respectively different from the classical 
partition functions found starting from cCHDa for similar chains. Thus, they differ 
in determinants of the sort referred to in Ref. [I]: QMa determinants are simpler 
than cCHDa ones. For ss macromolecules, we compare several quantities (bond-bond 
correlations, squared end-to-end distances, etc) from QMa with the standard Gaussian 
model in Polymer Science: the comparisons display good consistencies (which are also 
met with cCHDa). For double-stranded DNA (dsDNA) macromolecules, the Z^s from 
QMa have structures which generalize those obtained by other researchers, and enable 
to study thermal denaturation. 

There are many excellent general references on macromolecular chains (or polymers), 



1 



2 



1. Introduction 



from different standpoints. A (relatively small) set of them is: [2-15 

The contents of this tutorial review complements the subjects presented in the 
following contributions: [T| [T&j - [20] . 

Numerical computations for macromolecules are extremely difficult to carry out, 
because of the enormously large number of degrees of freedom involved. Fortunately, 
a subset of those degrees of freedom turn out to be, very approximately, constrained 
(taking on essentially constant values), namely, bond lengths, bond angles between 
successive bond vectors and, depending on the chain, other coordinates as well. The 
interest of accomplishing an adequate treatment of constraints in large molecular chains 
is easy to understand: it yields a reduction in the number of effective degrees of freedom 
to be treated and, hence, to a simplification in the simulations. The following pattern 
will be met in all macromolecules to be treated here. If the macromolecule is formed 
by N(^> 1) atomic constituents, one would have to deal, a priori, with aN degrees of 
freedom, a being some positive integer. By treating suitably various constraints, at the 
end, the number of effective degrees of freedom which remain unconstrained and, hence, 
have to be considered in the numerical computations will be bN, with b a positive integer 
such that b < a. 

The fact that the constituents of a macromolecule are atoms (described through 
Quantum Mechanics) raises some issues of principle, even if, at the end, one employs 
Classical Statistical Mechanics for practical purposes. Could one start the analysis of 
a constrained macromolecule safely from Classical (Statistical) Mechanics or, rather, 
should one begin with Quantum (Statistical) Mechanics and, at a later stage, proceed 
to its classical limit? The following argument will indicate that the answer to that 
question may not be straightforward. If, in a system of classical particles (say, strictly 
in the framework of Classical Mechanics), some coordinates (denoted by q CO n) are fixed 
(constrained), while the others (named as q) are not but evolve in time, then, as 
remarked by Brillouin 21 , the momenta canonically conjugate of q con will not vanish 
and will vary with time in general, in a way completely determined by both q and 
the momenta canonically conjugate to the latter. On the other hand, if, in a system 
of microscopic particles described trough Quantum Mechanics (for which, in principle, 
Brillouin's remark does not apply), certain coordinates q con take on essentially fixed 
values (with very small, or almost vanishing, uncertainty), then, according to the 
quantum-mechanical uncertainty principle, their canonically conjugate momenta take 
on any possible value and have an almost infinite uncertainty! 

We shall remind some works devoted to analyze constrained motions of simpler 
microscopic systems from first principles: the quantum-mechanical motion (via a 
Schrodinger equation) of a three-dimensional microscopic particle, which is constrained 
(by the action of some adequate potential) to move, in a suitable limit, along a given 
curve or on a surface 22-24 . The analysis was generalized in 25 to the quantum motion 
of a system described by a larger set of coordinates (qj and q CO n,i, to be unconstrained 
and constrained, respectively) [25]: under the condition that no crossed term containing 
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second partial derivatives (like [d/dqj] [d/dq con: i]) appears in the kinetic energy part 
of the quantum hamiltonian, then the quantum constrained system is described by a 
unique Schrodinger equation (which generalizes those for the surface and curve cases). 
Unfortunately, the kinetic energy operators typically associated to even the simplest 



macromolecular chains (with N > 3) will contain such crossed terms [see Eqs. (3.2) 



and (3.5)]. Then, the above condition is not fulfilled and there is no guarantee that the 
procedure in [25] would yield a unique quantum hamiltonian. See also comments in |26| . 

The above discussions contribute to confirm the inherent difficulties involved in the 
study of constraints in quantum systems. Is it possible to formulate a complete quantum 
theory of constraints for macromolecules which be free of difficulties? 

The following general aspects which will apply for all single-stranded (ss) 
macromolecular chains to be treated here. In three-dimensional space, we shall consider 
one single ss macromolecular chain formed by N non-relativistic atoms (treated either 
as classical or as quantum particles), with masses Mj(> 0). Let Ri be the position 
vector of the i-ih. atom, with i = 1, . . . , N, and let Mt t be the total mass of the whole 
ss macromolecule (see figure |i~T] ) . Also, let the center-of-mass (CM) position vector and 
the relative ones for the successive atoms (the bond vectors) be denoted by Rcm, and 
yi (where 1 < i < N — 1), respectively. One has: 



Rcm = ^P^ ,M tot = J2Mi, (1.1) 

i=i 

y 4 = R m -R. (1.2) 

Using spherical coordinates, we write: 

Yi = ViU-i , (1.3) 

Uj = (cos ifi sin#j, simpi sin#j, cos 9i) , (1-4) 

uq. = (cos ifi cos 9i, simpi cos 9{, — sin 0i) , (1-5) 

u w = (- sin<£i,eos</?j,0) . (1.6) 

Notice that, for fixed i, the vectors Uj, u ( p i ,ug i constitute an orthonormal set. 

The total number of atoms in any single macromolecule to be studied here (namely, 
N in the ss-case) is supposed to be enormously large. This justifies the applicability 
of Statistical Mechanics (and, eventually, of Thermodynamics) to those chains. The 
possibility of applying the latter frameworks to molecular chains which are not large 
has recently attracted attention (in particular, related to single- molecule stretching 
experiments) [27j[28]; we shall not analyze this subject here. 

We shall always deal with large three-dimensional macromolecules specifically in 
thermodynamical equilibrium, at absolute temperature T in an interval of physical 
interest, which includes room temperature (~ 300 K). We shall anticipate two ubiquitous 
physical facts. At T not far from 300 K, there are many allowed states for the 
macromolecule and, for each allowed state, all interatomic distances yi, i = 1, . . . , N — 1, 
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Figure 1.1: Single-stranded (ss) macromolecular chain. 

turn out to be approximately constant and ~ d% > (the bond lengths), namely, 
approximately constrained. Also, the energies associated to individual unconstrained 
angular degrees of freedom of microscopic constituents in the macromolecule are, 
typically, smaller than ksT (kg being Boltzmann' s constant). 

We shall restrict in Chapters [2] to [5] through to what characterizes the freely-jointed 
(or unhindered) chains. In them, all yi = di, i = 1,...,N — 1, while all angular 
variables will remain unconstrained (except for a closed-ring constraint in [5]) . Chapter 2 
outlines approaches to constraints in freely-jointed ss macromolecules, through Classical 
Hamiltonian Dynamics (cCHDa), displays certain controversies arising there and turn 
to quantum approaches for them. Chapter [3] justifies the quantum-mechanical physical 
assumptions and the variational approach to constrained macromolecules to be employed 
later. Chapter [4] treats freely-jointed macromolecules quantum-mechanically. Chapters 
[5] and [6] summarize generalizations for single-stranded macromolecules with constrained 
bond lengths and with further constraints (closed-ring and angular ones). Chapter [7] 
is devoted to double-stranded (ds) open chains with constraints. Chapter [8] deals with 
applications to dsDNA. 



Chapter 2 

Macromolecules with rigid and 
stiff constraints 



2.1 Classical models, classical partition functions and 
paradoxes 

Let all particles be regarded as classical. A dot above a time-dependent function will 
denote its first time derivative. The classical kinetic energy T a = | YliLi MiRf of the 



ss macromolecule, by employing Eqs. (1.1) and (1.2), also reads: 

N-l 



T - -M tot R 2 CM + 

*J=1 



\ J2 y l (B~ 1 hy J , (2.1) 



where the total mass Mtot was given in ( 1.1 ). Regarding the (N—l) x (N— 1) matrix B , 



with elements (-B _1 )jj, it suffices to know that: i) B~ l is symmetric {{B~ l )ij = (B ' 
ii) has positive eigenvalues (T a being a positive definite quadratic form); hi) the 
inverse matrix B, with elements Bij, also has positive eigenvalues and it is a tridiagonal 
symmetric matrix (for i = 2, . . . , N — 2, B^ = 0, unless i = j — + 1, and 

B13 = ••• = B\n-i = 0, -Bjv-ii = ••• = Bn-in~3 = 0). Explicit computations for 
low values of iV allow to confirm i)-hi) easily. See [29 and references therein about how 
to find analytically the inverse of a general tridiagonal matrix. 

We shall start with the constrained Classical Hamiltonian Dynamics approach 
(cCHDa), in which all yi = di{> 0), i = 1, . . . , N — 1, hold (holonomic or classical 
rigid constraints ). In so doing, there is not a unique procedure (alternative procedures 



being equivalent to one another). One, due to Kramers 30 and used in [31f|33] , employs 
only independent unconstrained variables. We shall disregard the overall motion of the 
CM. In Kramers' rigid model the independent unconstrained variables are all angles 
Qi and cpi, i = 1, . . . , N — 1, to be denoted collectively as qi, i = 1, . . . , 2(N — 1). 



As ji = di(6i\ig. + ipiUip.), with i = 1, . . . ,N — 1, Eq. (2.1) becomes (with the CM 

h0- 



disregarded): T 0jC(m = 2" 1 Y,fL l) G ijqi q r The [2{N - 1)] x [2(N - 1)] matrix G 
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with elements Gij, is obtained upon replacing the above y, into Eq. (2.1). Notice that 
Gij depends on qi, with i = 1, . . . , 2(N — 1). There may also be some potential energy 
V con = V CO n(qi, ■ ■ ■ , <?2(iV-l))> so t na * the lagrangian (with the CM disregarded) is L con = 
T 0yCon — Vcon- Throughout this work, the subindex "c" in a dynamical variable will always 
remind that the latter is classical (namely, not a quantum operator ). The classical 
momentum canonically conjugate to Qi and ipi are ne-i^ and 7r» ; i c , i = 1, . . . , N — 1, to 
be denoted, collectively, as 7r QrMjC . One has: n an ,i,c = dL con /dqi = Y?j=i Gijij- Tne 
associated classical hamiltonian is: 



H, 



2(JV-1) 



2(JV-1) 



(G- 



ij^an,i,c 11 an,],c 



7T„ 



(2.2) 



The matrix G _1 , with elements (G _1 )ij, is the inverse of G. The very large chain in 
thermodynamical equilibrium at temperature T is now described by Classical Statistical 
Mechanics, thereby disregarding quantum effects. We shall disregard the CM motion 
and concentrate on all qi and Tr a n,i,c- The classical equilibrium partition function for the 
macromolecule is: 



[du] [d7T a 



(27^)2(^-1. 



exp 



N-l 



2(N-l) 



[du] 



Y[ dOidifi , [d-n-an^] = ^| 



d-K„ 



(2.3) 
(2.4) 



By performing all Gaussian integrations on ir, 



an,i,cj 



Z C on = {"inksT) 



N-l 



[du] 



(27T^) 2 (^- 1 ) 



one finds: 



[det G] 1/2 exp 



V„ 



k B T 



(2.5) 



det G being the determinant of the matrix G. One generic difficulty of this classical 
model is that det G depends on the angles. See Fixman [34] for a study of det G. 
Other cCHDas to macromolecules ( equivalent to Kramers' one) work with all variables 
(all yi and their canonically conjugate momenta) and suitable Lagrange multipliers, 



see 126 , 35 . Another tutorial review in this volume |18| discusses generic iV-dimensional 



Lagrangian systems with D independent holonomic constraints, including the case in 
which the N — D unconstrained coordinates cannot be constructed explicitly, not even 
locally. The latter formulation (the ambient space one) and the associated Langevin 
and Fokker-Planck equations (to treat dynamics) are reviewed in [18| . 

The classical and unavoidable time variation of the momenta canonically conjugate to 
the j/j's (2l] turns out to be one of the great difficulties met in simulations employing the 
molecular dynamics method for the study of macromolecules, when constraints are taken 

Nevertheless, considerable progress along this research line has 



into account 136 37 



been facilitated mostly due to efficient algorithms allowing to integrate numerically the 
3iV Cartesian equations of a system of N classical point particles subject to holonomic 



constraints 38 . As a matter of fact, the holonomic constraints considered in [38] 
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included not only bond lengths but angular variables (bond angles) as well. The 
computational algorithms in [38 j (implementing molecular dynamics simulations and 
referred to, at present, as "the Shake ones") have been subsequently followed and 



developed further by a good number of authors. See 19 



By taking 35 as starting point and through a detailed analysis, Mazars [39 40 has 



obtained that the critical exponents for Kramers' model are the same as those of the 
Gaussian model [7 ( 13 . See also 41 . For other studies, aimed at comparing Kramers' 
model with the Gaussian one, see [26] and references therein. 

We now turn to another classical model considered by Fraenkel [42] (and pursued 
by |43}]53] ) , which is also based upon Classical Mechanics and employs harmonic springs 
(hs). Fraenkel's model, which is not equivalent to Kramers' one, is characterized by the 
fact that the constraints yi = di, with i = 1, . . . , N — 1 are not imposed from the very 
beginning, but at a later stage, by letting the harmonic springs to become very stiff: 
a classical flexible model with infinite stiffness (also known as stiff model). Flory has 
employed similar ideas in his classic work about statistics of macromolecules [4|. The 
lagrangian is now L^s = T — Vq, where T Q is given in Eq. (2.1 )) and the potential Vq now 
depends on all y\, 6i and <pi, i = 1, . . . , N—l. Upon introducing the canonically conjugate 
momenta 7Tj iC = dLhs/dji = Y^j=i (B~^)ijy j and disregarding the contribution of the 
CM, the actual hamiltonian reads: 



7V-1 



H, 



hsAn 



Hin + Vq, Hi n — - ^ 7Tj jCJ B 



(2.6) 



Bij are the elements of the matrix B [the inverse of B in (2.1 )]. The partition function 
for the chain at thermal equilibrium is: 

[dy] [dvr, 



Hhs. in 

W '* = ' (2^)3(^-1) 6XP [~1<bT_ 
JV-l AT-1 

^y] = n d3 y- t d7r J = n d37r ^ • 

i=l i=l 
Upon performing the Gaussian integrations over all momenta, one gets: 



(2.7) 
(2.8) 



Zhs 



(27rk B Tf( N - i y 2 



[det B] 



1/2 



[dy] 



(27r^) 3 ^- 1 ) 



exp 



Vr 







k B T 



(2.9) 



det B (the determinant of the matrix B) is a constant. We shall suppose that 



Vq = V con + U, where V con is similar to the one in Eq. (2.2) and in Eq. (2.5), while U 



is the following harmonic-oscillator-like (or harmonic spring) potential: 



N-l 



U = l J £u](B jj r\y j -d j ) 2 , 



(2.10) 
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where ujj(> 0) are frequencies and dj(> 0) are the bond lengths. The statistical average 
of a function F = F(y\, . . . , yjv-i) is: 



J[dy]F(yi, . . . ,y N -i)exp 


Vcon 

k B T_ 


J[dy] exp 


k B T 





(2.11) 



We shall assume that all Uj —> +oo (say, the limit in which the harmonic springs 
become very stiff), which forces exp[— [//(fc^T)] to equal a constant factor A" 1 times 
n^i 1 3{yj — dj), where 8 denotes Dirac 's delta function. Thus, the physical bond-length 
constraints yj = dj are recovered in the stiff harmonic spring limit. A should diverge as 
the frequencies do, although that will not be relevant here. Then, one gets in the stiff 
harmonic spring limit: 



dfi] exp 
dO] = Y\ sin6 s d6 s d<p. 



Vcon 



/[dft]Fexp 


Vcon 

k B T 


J[dft]exp 


k B T_ 



N-l 



(2.12) 
(2.13) 



s=l 



All constant factors in Zh s have been embodied into a single one, Ai, which, in turn, 
cancels out in (F). Those coordinates (like the y s ) which, as a result of the above stiff 
harmonic spring limit, take on constant values, are named hard variables [54] , The 
remaining coordinates, which remain unconstrained after having taken such a limit, are 
called soft variables [54]. The larger the frequencies are, the more rapidly and wildly 
oscillate the hard variables (y s ) about certain values (d s ), to be regarded as constant 
parameters. In any (small) period corresponding to those rapid oscillations, the soft 
variables do not change appreciably. Let V con = and let all d s = d, s = 1, . . . , N— 1, in 
Eqs. (2.12). In the stiff harmonic spring limit, Eqs. (2.12) yield: (y^) = d 2 , (u s Uj) = 
and ((Rjy - Ri) 2 ) = (N — l)d 2 . These results agree with well known results for the 
standard Gaussian model [5j [7j. A more detailed comparison shows that the three- 
dimensional distribution function for the end-to-end vector Rjy — Ri in Fraenkel's 
classical model with stiff harmonic springs coincides with the end-to-end distribution 
given in Eqs. (14-83), (14-84) and (14-86) in McQuarrie [6j. However, it does not seem 



that Eqs. (2.12) imply the Gaussian distribution for the individual bond lengths. See, in 
this connection, the discussion in pages 22 — 27 in (§]. Then, Fraenkel's classical model 
with stiff harmonic springs is partially consistent with the standard Gaussian model. 

We now turn to display paradoxes, upon comparing Kramers' model and Fraenkel's 
one. In fact, Z con ( (2.5)) and Zh s ((2.12)) differ clearly from each other. The discrepancy 
has been displayed neatly for suitable angular probability distributions in the case N = 3 
and V con 



(the trimer) [55-57 



In particular, the analysis in 55 was based on 
molecular dynamics simulations for N = 3. That discrepancy gave rise, in the framework 
of Classical Mechanics (without introducing, as yet, quantum effects and the uncertainty 
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principle!), to a controversy. The latter could be stated as follows. If, in the framework 
of Classical Mechanics, certain coordinates in the macromolecular chain ( the bond 
lengths) have to be regarded as constrained, should one employ Kramers' approach (rigid 
constraints) or Fraenkel's one (infinitely stiff flexible constraints)? Certain simulations 
favoured Kramers' approach, while others were more consistent with Fraenkel's one and, 



finally, others gave almost undistinguishable results. See 55 58 59 



2.2 On the role of quantum effects 

Several researchers carried out further analysis aimed at clarifying the above controversy, 
by bringing some quantum-mechanical ideas into the scene. Helfand [60] recognized 
clearly that the statistical properties of rigidly constrained systems (Kramers') and 
flexibly constrained ones with large stiffness (Fraenkel's) are different, in general. He 
also treated a simple two-dimensional quantum-mechanical model. See also Van Kampen 



57 . We shall focus below on 54, 61 and, mostly, on |56|. 



Go and Scheraga 54, 61 explored the origin of the differences between the classical 



partition functions resulting from Kramers' and Fraenkel's models [say, between Eq. 



(2.5) and Eq. (2.12)], by including from the outset quantum effects for the variables 
which should be constrained (the hard ones): in their analysis, the latter were not only 
the y s but also the bond angles. Their first work [54] supported Kramers' approach. 



However, the second work 61 concluded, that Fraenkel's approach was, within the 
approximations involved, more accurate than Kramers' one. 

Rallison [56] argued that the expectation that rigidity or rigid rods (Kramers ) and 
flexible stiffness (Fraenkel) be equivalent need be not valid in Statistical Mechanics, 
and that Quantum Mechanics is required to formulate the problem. According to him, 
the paradox had arisen from a failure to recognize the relevance of quantum effects, 
in the limit where certain coordinates are to be constrained to constant values. He 
also argued that a unique answer could be obtained in principle (although it could be 
very difficult to achieve, in practice), by regarding the system as the classical limit of a 
quantum- mechanical one. He carried out a detailed analysis for the trimer. According 
to him, the corresponding result for the trimer through Kramers' model would never 
appear as a natural limit, when a full quantum-mechanical analysis is performed for 
such a system with quantized vibrations. Accordingly, let oj s be frequencies, playing a 
role similar to those in Fraenkel's model. Rallison also tried to characterize the extent 
of quantum effects by the dimensionless parameter Q = hajg/ksT (K being Planck' s 
constant), which is reasonable, and stated that, when Q — > 0, the classical stiff spring 
results are recovered. The last statement leads to a conflict with the stiff harmonic 



spring limit in Fraenkel's purely classical framework, namely, when going from Eq. (2.9) 



to Eqs. (2.12) . Thus, coming back to the formulation of Fraenkel's model after the 
above comments, one now realizes that u s should not exceed ksT '/h, if the formulation 
is to remain in the classical framework, so that to allow for lv s 3> ksT/h would bring 
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in quantum effects. In [56], the degrees of freedom to be constrained (due to the very 
stiff springs) were treated quantum mechanically, while the remaining unconstrained 
coordinates were regarded as classical. 

No strictly definite conclusion about which model was more accurate seemed to 
have been reached. Kramers' model is formulated in a mathematically consistent 
way in the frameworks of (constrained) Classical Mechanics and Classical Statistical 
Mechanics. A number of researchers have made the following choice: they have 
disregarded Fraenkel's model, and they have concentrated on developing quite vigorously 
computational approaches based upon constrained Classical Hamiltonian Dynamics 



(cCHDa): see |37[ 38 and other articles in this volume. On the other hand, other 



researchers tended to adhere, with various qualifications and reservations (mostly arising 



from quantum considerations), to the conclusion in 61 , namely, Fraenkel's approach: 
see Section 1 of |58| (where Brownian dynamics simulations for classical chains are 
used). For a comprehensive discussion, see the Introduction in [62], which includes 
a lot of additional and updated information. The different effects of the constraints 
in both Kramers' approach and Fraenkel's one, on the conformational equilibrium 



distribution for a relevant biomolecule (a dipeptide), were studied comparatively in 62 
although the question of which model was a better approximation was not specifically 
addressed. This work analyzed for the first time the conformational dependences of 
correction terms (related to the determinants) to that equilibrium distribution, using ab 
initio quantum mechanical calculations (including electron correlations) and without 
simplifying assumptions. Their numerical computations (Monte Carlo simulations), 
which employed a realistic potential energy function, concluded that those correction 
terms could be neglected in certain situations, up to peptides of considerable length, 
but not in general. They also analyzed whether the equilibrium values of the hard 
coordinates were dependent on the soft ones: their computations indicated that such a 
possibility should be critically considered in each case. See [l]. 

What about quantizing Kramers' model? The classical Kramers' model for a 



constrained freely-jointed molecular chain, based upon Eq. (2.2), can be quantized 



through certain well established procedure initiated by Podolsky 63 and developed 



further by De Witt 64 , with a probabilistic interpretation in the standard quantum- 



mechanical framework [65, 66 . For details, see 26, 67 and references therein. Since the 
quantization through |63| is formulated, from the outset, in terms of the unconstrained 
coordinates (say, 9 S and ip s ), the former will shed no light on the issue of the quantum 
uncertainties for the constrained coordinates (y s ) and of those for their canonically 
conjugate momenta! 

One could also follow the alternative path in which Fraenkel's model is not 
disregarded and to proceed further, by entering definitely into the quantum regime. 
In such an option, it is worthwhile to extract the following summary from the comments 



and studies above. Thus, in agreement with [54 56 61 , we remind that: i) the best 



description of the actual molecular chain is provided, in principle, through a quantum- 
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mechanical treatment, ii) specifically, one should start with a quantized chain with 
flexible but stiff springs, iii) at a later stage, one could proceed to the classical limit. We 



also recall that the main troubles regarding this philosophy were practical 54 56, 61 



it was very difficult to carry through a treatment based upon i) and ii) , unless a number 
of approximations were made, the validity of which, in turn, was rather difficult to 



control. Thus, as commented in 54, 56, 61]: a) the large quantum zero-point energies 
of the constrained hard degrees of freedom (due to the very stiff springs) could depend, 
in principle, on the remaining unconstrained (soft) coordinates, regarded as classical, 
b) if such dependences were taken into account, the subsequent analysis appeared to 
be quite hard, while if they were not, then the validity of the resulting approximations 
would be difficult to assess, c) if those quantum zero-point energies were sufficiently 
large (as they would eventually become, for suitably large frequencies of the harmonic 
springs), it would be necessary to treat, at least in principle, the remaining unconstrained 
coordinates not just as classical variables, but through Quantum Mechanics as well: the 
latter possibility could add even greater difficulties. 
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Chapter 3 

A variational 

quantum- mechanical approach: 
General aspects 



What about quantizing Fraenkel's model? If done adequately, do large quantum zero- 
point energies of hard degrees of freedom depend on soft ones? 

The models for macromolecular chains based on Quantum Mechanics to be reported 
in this work will be derived through a different (and less ambitious) strategy: the 
application of a variational quantum-mechanical inequality procedure. The latter, not 
relying on the absence of the crossed terms containing second partial derivatives, will 
not be subject to the limitation in [25]. Their variational foundations will imply that 
the quantum models reported here should be not be regarded, in principle, as definitive 
or final formulations, but as approximate ones. 

The following comments should suffice to support physically the variational approach 
to be pursued along this paper. The physical fact that the relative distances 
(bond lengths) from any atom to its nearest neighbours are, within narrow limits, 
approximately constant, in generic macromolecular chains (below some maximum 
temperature ) I2H7, 11 , provides justification to the following picture, which incorporates 



the Born-Oppenheimer approximation 68 and the previous remarks |54l 56, 61 



All atomic constituents in the macromolecule are subject to an effective potential 
(due, reciprocally and self-consistently, to all of them), which includes: 1) electronic 
contributions (the largest ones, denoted as E e i), 2) potentials (having magnitudes smaller 
than those in 1)) which, for intermediate distances, resemble qualitatively harmonic- 
oscillator-like vibrational ones and force atoms to oscillate with frequencies Ui about 
some equilibrium position in the chain, so that the relative separations to its nearest 
neighbors do not vary appreciably, on the average, and 3) other potentials (having 
magnitude smaller than the vibrational ones and, hence, not upsetting the oscillations) 
and are, in turn, responsible for other weaker effects: angular variables, etc. As another 
justification, we recall that for certain simple (small) molecules about room temperature, 
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individual internal rotations about bonds are slowly-varying, as their typical energies 
are, at least, two orders of magnitude smaller than the vibrational energies hui [3 68 



such a qualitative pattern appears to be valid for macromolecules as well, on the average 
(except, possibly, for some moderate quantitative variations). For generic vibrational 
frequencies u;, of typical atomic constituents in physically interesting macromolecular 
chains, it is a fact [3] that the vibrational energy huii is larger (about one order of 
magnitude more or less) than ksT (T in an interval about 300 K). Thus, for a chain at 
thermal equilibrium at those ksT, generic vibrational degrees of freedom are rapidly- 
varying, have to be treated quantum-mechanically (in a way radically different from 
the classical description), and they appear to be, quite approximately, in their ground 



states, with all excited vibrational states unoccupied 54, p6|. Electronic degrees of 



freedom remain unexcited in phenomena involving energies in the above range of /cgT 
and, hence, E e i can be regarded as approximately constant. 

The above facts agree with certain qualitative statements made by Schrodinger in his 
well known book |69|: in fact, his statements (on quantum theory, molecules and how 
their stability depends on temperature) appear to have discussed very lucidly a good part 
of these physical issues. It follows that quantum effects for those vibrational degrees of 
freedom will be dominant (and even more for electronic ones). Hence, one should employ 
Quantum Mechanics as the starting point to deal with them, at least in principle. Those 
dominances and the fact that electronic and vibrational energies are larger than 
(T in an interval about 300 K) are ubiquitous manifestations of quantum effects in 
macromolecules (in particular in biological ones). In turn, the quantized vibrations 
about their ground states in macromolecular chains at thermal equilibrium appear to 
decouple, at least as a zeroth approximation, from internal rotations. Our study, with 
uj s S> ksT/h could be viewed as a fully quantized version of Fraenkel's one. It will also 
improve that in I56J , since all variables (those to be constrained and the ones which will 
remain unconstrained) with be treated quantum-mechanically from the outset. 

The variational quantum-mechanical models to be treated from section [4] onwards 
would be subject, in principle, to the objections a), b) and c) at the end of subsection 



2.2 We anticipate that, fortunately, by virtue of various crucial exact cancellations, at 



the end of the computations those models will turn out to be free of the difficulties a), 
b) and c) at the end of subsection 2.2 the expected quantum zero-point energies of hard 
degrees of freedom will be obtained and, hence, shown to be constant. 

The successive analysis of quantum macromolecules with constrained bond lengths 
and increasingly complicated constraints are neither easy nor direct. A number of those 
generalizations have indeed been carried out, but at the expense of getting involved 
into and overcoming considerable difficulties at the quantum level. One finds new 
and interesting quantum-mechanical structures. Although it is not easy to handle 
those structures, they imply an important consequence: the rotational invariance ( 
conservation of total angular momentum) in the quantum models. Moreover, those 
structures provide a consistent basis to get classical approximations, at about room 
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temperature, which are more manageable. 

We shall employ both harmonic-oscillator-like potentials and Morse potentials 



[66 70 . The latter can be approximated by ( and are less idealized than) the former, 
for certain range of intermediate distances. We shall outline the methods and avoid 
the presentation of calculation details throughout the main body of the present tutorial 
review. Some specific calculations are given in the Appendixes. 

We will devote a great deal of attention to constructing partition functions for 
the different chains, which will be well defined and meaningful first at the quantum 
level and, later, at the classical one. We shall carry out several approximations in 
the more relevant parts of the resulting classical partition functions. Specifically, we 
shall use those classical partition functions (with or without those approximations) in 
order to find various physically interesting quantities: internal energies, correlations, 
probability distributions, equilibrium constants. For the latter purposes, the less relevant 
contributions to those partition functions will not contribute and, hence, will not be 
given. Thus, the complete results of the approximate evaluations of the classical 
partition functions will not be strictly needed and, hence, will not be reported: see 



26 71-75 



3.1 The variational inequality 



We continue to use Eqs. 
operator: 



( |1.1[ )-( JL6"| ). We start from the total quantum Hamiltonian 



H, 



re 

2 



N 



£ M < + V(R), 



(3.1) 



where V denotes the gradient operator and V(R) represents the real potential energy 



among atoms. The kinetic energy operator — (^ 2 /2)^^ 1 M i 1 Vj l . in Eq. (3.1) is the 



quantum- mechanical version of the classical kinetic energy in Eq. (2.1). Then, using 



expressions (1.1) and (1.2), (3.1) gives the quantum Hamiltonian in terms of the new 
variables as H Qtl = -{h 2 /2M tot )Vl + Eq. CM de grees of freedom will always be 
factored out and disregarded (and, hence, so will be — {h 2 /2M to t)^ 2 \ CM ), for all types 
of single-stranded macromolecular chains to be studied in Sections [3](6} But CM degrees 
of freedom will be taken into account for double-stranded chains (Sections [7] and [8j. For 
single-stranded chains, we shall always concentrate on the internal quantum Hamiltonian 
operator: 



H. 



Q 



2 



N-l 



N-l 



Mi 



+ U(y) = H Q>in + U(y) 



(3.2) 



Hq^u is the quantum-mechanical operator associated to 
coefficients Bij in Eq. (2.6) are those in HQ^ n : in particular, Bi 



in Eq. ( j2J3| ). The 

= M- 1 + M-\(= J Ey. 
See also Eq. (|4.26|). For single-stranded chains, we denote the set of atomic coordinates 
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(yi) y2> • • ■ 3 yjv-i) by y (for double-stranded chains, in Sections 7JJ8, y will have another 
meaning). We remark that all molecular chains are treated in the framework of the 

Then, by recalling the previous discussion 



Born-Oppenheimer approximation |68, 71 



in this section |3| the effective potential in the macromolecule is described by E e t + U(y). 
E e i(< 0) is the electronic energy, corresponding to the electronic degrees of freedom (the 
most rapidly- varying and the most energetic ones), which is regarded, essentially, as a 
constant. We shall suppose that E e i has already been subtracted out and, then, omitted 
from the outset. There remains a total (real) effective potential energy U(y) = V(R), 
independent on Re m, which accounts for the effective interactions denoted as 2) + 
3), previously in this section. U(y) accounts for large covalent-bond interactions, 
which are responsible for the very existence of the macromolecule, as an extended and 
connected object, with energies smaller than the electronic ones (E e i), and other weaker 
interactions. Specifically, U(y) includes, first of all, the interactions which constraint any 
I yi |= y« to equal the fixed bond lengths, as well as the (somewhat weaker) interactions 
hindering rotations (constraining bond angles), those associated to the unconstrained 
angular variables and further weaker residual interactions. 

The three-dimensional momentum operator reads in spherical coordinates: 

d 



iHu, 

Vi oyi 



i 



d 



(3.3) 
(3.4) 



and we recall Eqs. (1.3)(1.4), (1.5) and (1.6). We shall refer to both yi and u n as bond 
vectors. For later use, we denote the set of all 6\, . . . , 9n-i, ¥>i> ■ ■ ■ > fN-i by 9 and ip, 
respectively. Notice that the a$'s are not Hermitian operators. Using Eq. (3.3), the 
atomic kinetic energy operator Hq^ becomes: 



Hr 



1 



N-l 

£2% 



1=1 

N-l 



-ft 



dyf Vi dyt yf 

|2 



V- 1 f,2 d 2 . fc 1 d \ 

ti. Mi { dyi-idyi yiOyi-i) 



i=2 
N-l 



Elf., Id a*_i • an 
— < ih&i-x ■ Ui 
M < { Vi-i oyi yi-\y% 



i=2 



(3.5) 



Since, by assumption, the large quantum (Q) macromolecular chain is in 
thermodynamical equilibrium, at absolute temperature T, we describe its state through 
the (exact) quantum partition function Zq, which is given by 



76 



Z Q = TrfexpH^T)- 1 ^]] = ^exp[--^ r } , 



k B T' 



(3.6) 



where a denotes the set of all quantum numbers and E a represents the eigenvalues of 
the entire spectrum of the quantum Hamiltonian Hq. 



3.1 The variational inequality 
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Let &\ be an arbitrary ortho-normal set of wave functions for the system (A 
representing a set of quantum numbers). It is not required that $a coincide with the 
exact eigenfunctions of Hq. Then, the exact quantum partition function for the system 



Zq (as given in Eq. (3.6)) satisfies Peierls' variational inequality 
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Z Q > Y,^P[-(k B T)- 1 (^ x ,H Q ^ x )] . (3.7) 

A 

The equality holds if <f>\ is the complete set of exact eigenfunctions of Hq. The <&\S will 
be the variational trial wave functions. The partition function Zq and the inequality Eq. 



(3.7) will provide the basic framework for the successive variational computations to be 
summarized here. We anticipate that, before performing the variational computation in 
each case, a very crucial and delicate question will be the choice of suitable trial wave 
functions (say, the $a's). For that purpose, the following key fact, genuine of the Born- 
Oppenheimer approximation |68| will be very helpful. On physical grounds, the (fast) 
vibrational and the (slow) rotational degrees of freedom may be decoupled, at least as 
a zeroth order of approximation, because the vibrational energies are larger than the 
rotational ones [68] . 

At about room temperature, rotational configurations of groups of atoms in ss 
macromolecules can be regarded to be, typically, in excited rotational states. This 
fact suggests that, at about 300 K, the slow rotational degrees of freedom could be 
treated through Classical Statistical Mechanics. Based upon all that, the successive 
approximations presented in sections 4, 5 and 6 for ss chains will follow the following 
pattern: we shall start out from a quantum-mechanical formulation, deal firstly with 
and decouple the quantized vibrational degrees of freedom, turn to the rotational ones 
(at the quantum level first) and take the classical limit thereof. At the end, the classical 
statistical description will appear as a justified approximation. Our final purpose will be 
to arrive, following such a route, at a simpler and reliable effective classical hamiltonian 
and partition function for ss macromolecules, depending only on the relevant slowly- 
varying degrees of freedom (all angular variables, in this case). That pattern will be 
extended to double-stranded chains in [7J 
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Chapter 4 

Single-stranded open 
macromolecules 



4.1 Freely-jointed chains: radial variational computation 



We shall treat in |4.1| and |4.2| freely-jointed chains (only the bond lengths being 
constrained). Freely-jointed chains constitute an idealization because, in addition to 
bond length, angles are also constrained in real macromolecules. In spite of that, the 
study of freely-jointed chains is absolutely essential, because the interactions which 
produce the bond length constraint are the most important contributions to U(y) and, 
hence, are responsible for the existence of the chain as an extended connected object. 

We assume that U (y) includes interactions among the atoms only if they lie in nearest 
neighbor positions. We shall consider two possibilities. In the first, U(y) describes 
harmonic-oscillator-like potentials (and equals U, for Fraenkel's model): 

JV-l 2 
i=l 

0Ji(> 0) are frequencies, while dj(> 0) are bond lengths (the equilibrium distances 
between two successive atoms). In the second possibility, we use the Morse potential |70|: 

V M (y) = D{exp[-2a(y-d)]-2exp[-a(y-d)}} . (4.2) 

D is the dissociation energy corresponding to the equilibrium distance d and a -1 



represents the range of the potential. Then, instead of (4.1), we consider: 

iV-l 



U(y) = E V MiVi) > (4-3) 



i=l 



each V~M{yi) with its corresponding a\ and D\. We shall introduce: uji = ai\j2DiBi. 



Our purpose is to arrive, through a variational computation [via Eq. (3.7)], at a 



model for a microscopic chain, in which all j/j be constrained at the quantum-mechanical 
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level. We shall follow the presentation in |72|, to which we refer for further details. By 
recalling the difference between fast vibrational (hard) and slow rotational (soft) degrees 
of freedom, we choose the trial variational wave function 3>(y) as a product of two wave 
functions <p{y) and tpa(0, <p) depending, respectively, on the radial and angular variables 
(cr denotes a set of quantum numbers). Our choice of the trial radial wave function is: 



JV-l 



<t>{y) = ]~[ Mvi 



(4.4) 



i=i 



Then, the (normalized) trial variational (radial plus angular) wave function of the chain 
reads: 



$ = $(y) 



"JV-l 



n ^ m 



i=i 



(4.5) 



Throughout this work, we shall suppose that all vibrational frequencies ui are very large 
(and, more or less, of similar orders of magnitude) and fulfill: 



huji > k B T, fvjji > d, 2 h 2 Bi 



(4.6) 



with the corresponding uj\ for either (4.1) or (4.3). For Eq. (4.1), each 4>i(yi) is chosen 
as a Gaussian 65 , approximating the ground state of the corresponding harmonic- 
oscillator-like potential: 

Uvi) = di l (t^w) exp \-T^B-(yi ~ dif 



I = 1,...,N- 1 . 



(4.7) 



In order to perform calculations to leading order in frequencies, we shall make use of 
the following property of Dirac's delta function (5(y — d)): 



lim 

A— »oo \ 7T 



1/2 exp [_ A(y _ d?] 
d 2 



S(y-d) 
d 2 



(4- 



66 



70 



For Eq. (4.2), we consider the normalized radial wavefunction <zW,n=o(y) 
associated to its ground state. It corresponds to the discrete ground state energy Em,ti=q'- 

hw h 2 u 2 i 

E M ,o = -D + — - ——, lo = a^/2DB J , (4.9) 



and so on for (4.3), with ground state energy YldLi 1 ^M,l,n=o- ^M,i,n=o is given by 
(4.9), with Di, ai, ui\. We shall notice the following useful property 

S(y - d) 



\<l>M,n=Q{y)\' 



d? 



(4.10) 



when a —> +oo and D — > +oo, while the dimensionless ratio y/2D / 1 {ha[Bj\ 1 / 2 ) equals 
a finite constant, so that uji = a\\j2DiB\ grows (say, when the Morse potential is very 
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deep). Eq. (4.10) plays for the Morse potential the same role as Eq. (4.8) for the 



harmonic-oscillator-like one. We have obtained a mathematical justification of (4.10) 



which is a bit lengthy and, so, will be omitted. An important feature of the Morse 
potential is that it is a bounded potential, so that the Schroedinger equation containing 
it has a finite number of bound states and it can account for both bound as well as 
unbound states. This clearly distinguishes it from the harmonic oscillator-like potential, 
which rise to an infinite number of bound states, but no unbound ones. 

Hence, as each ui — > oo (I = 1, . . . , N — 1), both Eqs. (4.8) and (4.10) yield: 



2 , TT S (Vl ~ d l) 



n 



(4.11) 



The normalization condition J[dy] | 3>(y) | 2 = 1 ( with [dy] given in (2.8)) becomes in 
the limit ui — > oo: 



[d«]|^(0,V>)| : 



1 



(4.12) 



with [dfi] given in Eq. (2.13). We emphasize that, except for the above normalization 
condition, Vv($ ; V 7 ) is fully arbitrary. We evaluate, as all frequencies Ui, i = 1, . . . , N— 1, 
grow very large, the quantum expectation value (<&,Hq&) = J[dy]$(y)*iJg$(y), by 



using Eqs. (4.11), (4.5) and (4.12). The details of the computation are summarized in 



appendix A. The result can be cast as: 



(VvX o) Vv) 



(4.13) 
(4.14) 



o and J2^ = i l (the sum of all 



for any ip a (8,ip). Here, Eq stands for Ya=i^ ^M,l,n 
zero-point energies) for all Morse and harmonic-oscillator-like potentials, respectively. 
The new quantum (angular) Hamiltonian Hq^ reads: 



H. 



(o) 



N-l 

E 



B; 



d? 



N-l 

E 



M; 



(ej-i ■ ej) 
didi-i 



(4.15) 



with e; being the following operator: e; = ihui — a^, I = 1, . . . , N — 1. Through direct 
partial integrations, e/ can be shown to be a Hermitean operator. Notice that Hq^ 
is expressed in terms of the e^'s, which could be regarded as some sort of quantized 
transverse momenta. Finally, 0^°'{K) denotes the set of all remaining contributions 
which do not depend on the frequencies oji : they are proportional to some positive power 
of h. Contributions of this sort will also arise in the analysis of other quantum chains, 
with more complicated constraints. In this tutorial review, we shall not give explicitly 
the various 0(H)'s which arise successively, because they are increasingly complicated, 



can be obtained from 26 71 - 73 75 and become eventually irrelevant when one takes 
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h — > 0. We shall concentrate on the Hq^s, due to their interesting structures and 
because they give a nonvanishing classical limit. 

It is important to notice that all dependences of (<3?, Hq&) on the frequencies appear 
only in the constant Eq, in the right hand side of Eq. ( |4.13| ). This solves the difficulties 



a), b) and c) at the end of subsection 2.2 



We shall state very briefly the following mathematical properties of the e^'s, which 
correct and complement adequately a previous discussion in [26], and which will be 
physically relevant at the end. Let us write e/ = (eji, e^, ^3), ej i and so on being the 
three Cartesian components of ej. Let 1; = y/ x [— iKV yi ] = (Zj i, k t 2,h,3) be the quantized 
Hermitean) orbital angular momentum associated to y;. Let [A, B] = AB — BA be 
the commutator of the operators A and B. The commutation relations of In, l\2 



and k t 3 are well known in Quantum Mechanics 65, p8j: [h,aJj,/3] = i^k,j^ai3^j,-y 



a, (3 = 1,2,3 65|68j. Here, S^j and e Qj g 7 are, respectively, the Kronecker delta and 
the standard totally antisymmetric tensor with three indices (£123 = +1, etc.) One 
could ask what are the commutation relations of the operators en, e\% and e^. Such 
a question is quite natural, at least from a purely mathematical standpoint, as an 
attempt to characterize the e^'s somewhat more. Through some lengthy algebra, one 
gets (a, P = 1, 2, 3): 

[efc,a> e,-,/?] = — ^fc,i e a/37^',7 > [4,m e j,/3] = ^fcj e a^7 e j,7 > (4-16) 

to which the above commutation relations of In, I12 and ^3 should be added. A 
simple, but important, consequence is that Hq^ is a Hermitean operator. One also 



has: e| = l| + ft 2 and e^.lj = Ij.ej = 0. Then, Eqs. (4.15), (4.16) and the commutation 
relations of Zj i, Z; ; 2 and Z/,3 imply, consistently, that the quantum total orbital angular 
momentum Ylf^ih commutes with Hq\ which is physically important. See 



2G 



and references therein. Thus, the conservation of the total orbital angular momentum 
fully justifies, a posteriori, the previous (painstaking) excursion into those commutation 
relations and algebraic matters. We shall add here the following curious feature (not 
commented before) of the above algebra for the quantum freely-jointed chain. The 



commutation relations in (4.16) together with those of In, I12 and li 3 coincide with 



the ones for the generators of the homogeneous Lorentz group, which, in turn, play a 



crucial role in Special Relativity 78 . Non-relativistic quantum-mechanical approaches 
to macromolecular chains and Special Relativity appear to be subjects very disjoint 
and disconnected from each other. So, the unexpected appearance of the same algebraic 
structures (the same closed algebra, formed by the above commutation relations) in both 
subjects is certainly curious and surprising. Then, Hq ^ and the above algebraic structure 
at each bond vector describe an extended, connected and flexible quantum system (some 
sort of non-relativistic quantum "string", loosely speaking). For the quantized version 
of Kramers' model, one can introduce a set of variables which, even if different from 
efc.o's, play a role similar to the latter. Such variables in Kramers' model, together with 
an adequate formulation of the corresponding quantum angular momentum variables, 
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can be shown to generate a closed algebra identical to the above one for efc jQ 's and k } ps. 



We omit details and refer to 67 



We proceed to implement Eq. 



3.7) in this case, with $ A = $(y) and (<&\,Hq$\) 



being given in Eqs. (4.5) and Eq. (4.13). Then, Eq. (3.7) becomes: 
Z Q > exp [-{k B T)- x E G ] ■ Z Q o) , 



Z, 



(o) 



J>xp -{kBTY 1 j[dilWA0^)H { Q ] M0^) + O^{h) 



(4.17) 
(4.18) 



Zq can be regarded as the approximate three-dimensional quantum partition function 
for the slowly- varying angular degrees of freedom of the open (o) chain (in the framework 
of Peierls' variational inequality). The bond lengths are constants and, so, constrained 
in the quantum-mechanical Zq . Notice that if O^(h) is disregarded and the angular 
wave functions Vv (fl, ip) are taken as the complete set of all orthonormal eigenfunctions 



of H Q °\ then Eq. (418) becomes Z Q °> = Tr[exp [-(k B T)- l H Q L 



(o) 



- X H%\ This Z { ° ] includes the 



J Q 



lj) 2 ) and of its third component (X^i 1 h 



contributions of all possible quantum states and, hence, for all possible values of the 
squared total angular momentum (Q^ =1 ± 3 
Those values being the standard discretized values 
are interesting researches on classical partition functions of single-stranded chains with 



3 =i -j,3J- 

65[ |66|. In recent years, there 



configurations corresponding to a given value of the total angular momentum. 79, 80 



The construction, through the above variational procedure, of quantum partition 
functions including only states with prescribed values of (%2f=xlj) 2 an d Tl!j=i h& 
stands as an open problem. 

Let all atoms be identical to one another (Mj = M in the chain, i = 1, . . . , N) and, 
for simplicity, let all bond lengths be equal ( dj = d, j = 1, . . . , N — 1). Then, the very 
existence of the bond length constraints, forcing identical atoms to lie successively along 
the chain destroys their indistinguishability to a very large extent, but not completely. 
A discussion for the quantum Kramers' model (for both bosons and fermions) appears in 
|67| , which appears to be also valid for the actu al (v ariational) quantization of Fraenkel's 
model. The trial variational wave function 



4.5) and its radial part ( n^Li* filivi) 



also hold, while ifj a {9,tp) has to fulfill the surviving indistinguishability restrictions. 



We shall apply briefly the argument leading to Eq. (4.18) to a one-dimensional 



quantum macromolecular chain, at thermal equilibrium. By assumption, the quantum 



Hamiltonian is given by the right-hand-side (rhs) of Eq. (3.2) with y, and V yi replaced 
by yi and d/dyi, respectively, the one-dimensional relative coordinate j/j being positive 
or negative. Also, the potential energy among the atoms in the chain is provided by 
the rhs of Eq. (4.1), with yi substituted by | ?/i |. The variational trial function is 
given by the rhs of Eq. (4.7), with d j 1 re placed by 2" 1 and y\ substituted by | y\ \. 
Then, the actual counterpart of Eq. (4.13) becomes, simply: (<&,Hq&) = Eq + 0(h), 
with Eq = YliL'i 1 ^ an d O(h) being frequency-independent and of order K, at least. 



Eq. (4.17) also holds, now with Z$ = exp[— {k^T) 1 0(h)} (which, in turn, approaches 
unity as h — > 0). Notice that exp[— {ksT)~ l Y^i^ ^ s partition function for a 
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set of harmonic oscillators in the high-frequency limit, and, so, it differs drastically from 



that for classical oscillators 65 



So, for large vibrational frequencies, we have obtained Hq from Hq in (3.2 ) and Zq^ 
from (3.6), through (3.7). We remind here certain mathematically rigorous results (as 
equalities in asymptotic limits, not as variational inequalities) for other different models 
by several authors: see 81-85 and references therein. Their results bear qualitative 
similarities to ours, although there are also considerable differences. In fact, those 
authors start from certain hamiltonians (depending on both hard and soft variables) 
and, without using the variational inequality ( |3.7[ ), take the large coupling constant 
limit (akin to our large vibrational frequencies), so as to obtain: i) asymptotic behaviors 
of (low- lying) eigenvalues of Schrodinger operators [81] , ii) simpler hamiltonians and 
partition functions for spin chains (Haldane-Shastry chains), spins being here the soft 



variables 82-85 



Let two external stretching forces f and — f act upon the two atoms at Ri and Rjv, 
respectively. Then, Hq in (3.2 ) is replaced by Hq? = Hq + f X^I 1 Yi- The quantum- 



mechanical variational computation goes through as above and yields a new quantum 
partition function: Zq\ = Tr[exp[— (ksT)^ 1 (Hq^ + f YliLl 1 <^ u i)]]- 



4.2 Freely-jointed chains: classical partition function 



Eq. ( |4.15[ ) and Eq. ( |4.18| ) simplify enormously, if one proceeds to the classical (h — > 0) 
limit for the angular degrees of freedom, under the additional assumptions: 



k B T > d7 2 h 2 B, 



(4.19) 



so that excited rotational states are occupied and quantum operators and statistics can 
be approximated by classical ones. Then, (h) and all quantities of order h or higher 
( and not c onta ining any uji) can be neglected, so that Eq. (4.6) be respected. Eqs. 
(4.15) and (3.4) indicate that the Hamiltonian Hq^ becomes, in the classical limit: 

N-l , 

b a ' 

8=1 



2d 2 



N-l , 

El \ a i,c ' a i— 1, 
M 

i=2 1 



did. 



(4.20) 



The aj !C 's are no longer operators but classical variables ( the classical limit of e^). They 
are: 



sin 9; 



(4.21) 



-Pfli) P<pi are the classical momenta canonically conjugate to 9i and (fi, i = 1, . . . , N — 1. 

The three-dimensional classical partition function (the classical limit of Eq. (4.18)) 
is (\dP 6 dP 9 \ = W^ismO^dPe.dP^) 
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Z (o) 



(27Th) 2 ( N - 1 ) 



dft] / [dPedP^ 



H 



(°) 



exp 



k B T 



(4.22) 



4.2 Freely-jointed chains: classical partition function 



25 



Are there further classical constraints ( besides y\ = di), contained somehow in 
(4.20)) and (4.22)), which would be the counterparts of those expressing iji = in 



Kramers' classical model?. Let us introduce, for a short while, the classical momenta 
Ttj^F = 7r j,c,F' a j ~ dJ 1 &j tC ( not to be confused with 7Ti )C , employed and integrated over 



in Fraenkel's classical model in 2.1 7Tj tCi p being a radial momentum. One can develop a 
mathematical argument implying that TTj lC ,F 



86 



The integrations over all the classical momenta in zjf^ are Gaussian. Upon 



performing them 86 , one finds: 



N-l 



(detB) 3 / 2 

1-1/2 



2nh 2 _ 

z R = /[dn][Ajv_ij 

[A JV _ 1 ]- 1 /2 = [detCuiCB- 1 )^-^-)]- 173 , 
The elements Bij of tridiagonal matrix B, of order (N — 1) x (N — 1), are: 



(4.23) 

(4.24) 
(4.25) 



B 



M, ^ M i+1 
1 

Mi 



if i = j 

if j = % — 1 or j = i + 1 
otherwise 



(4.26) 



On the other hand, the matrix (ui(B 1 )ijUj) is given by: 





r i 






E^ + i 




if i = j 




1 

Mtot 


ELi Af* 




E£ i+ i M 


( u r ; • U rj ) 


if i < j 




1 

[ Mtot 


Ei=i W 






(u r , • U rj ) 


if i > j 



Mtot = EaLi A^fc and (ui(-B -1 )i 3 -Uj) is a symmetric matrix of order (N — 1) x (iV — 1). 
The determinant [A^v-i] -1 / 2 is manifestly rotationally invariant (as it depends on the 
scalar products (uj-Uj) but not directly on angles). However, the matrix (iij(.B )jjU,) is 
not tridiagonal, unlike (Bij). The contribution of [Atv-i] 1 ^ 2 differs from those of other 
determinants appearing in [TJ |30j[34j [39} [61} [62] . The equilibrium statistical average of 
a function F = F(\ii, Uj\r-i) is: 



(4.27) 



Z, 



(<->) 



2vr/i 2 



JV-l 



(detB) 3 / 2 



dm 



F(u) 



[A^_x]V2 



(4.28) 



The internal energy U of the linear polymer, which is an interesting thermodynamical 
property, does not require the evaluation of [Ajv^iJ^ 1 / 2 . Thus, U is obtained from 
U = A — T(dA/dT), where A is the free energy (Z^ = exp[— (ksT)~ l A]). One finds 
readily energy equipartition, namely: U = (N — l)ksT [26} [72}. 

We shall suppose that the region in which the ss open freely-jointed chain moves 
(inside which its center-of-mass Rcm varies) is a sphere of very large radius Rq. Let 
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Mi = M (with i = l,..., N) and dj = d (with j = 1, . . . , N — 1), namely, equal masses 
for all atoms and equal bond lengths in the chain. Let us introduce the probability 
distribution function for the position vectors Ri,. . . ,Rat: 

W(Ri, ...,IL N ) = (T[ (R l+1 - R, - du^ , (4.29) 

5^ being the three-dimensional delta function. Similarly, the probability distribution 
function for the end-to-end position vector is: 



N 



R 



<S< 3 > R 



V 



R 



N-l S 

i=i / 



(4.30) 



We remind the standard Gaussian model for ss open freely-jointed chains: 

N-l 



^ e9 (Ri,...,Rjv) 
W G (K l+1 -K r ,2d 2 ) 



eq 



Y[ W G (-R l+1 --R f ,2d 2 



i=i 

72\l3/2 



exp[ 



R 



/(2d 2 )] 



(4.31) 
(4.32) 



Wg(Ri+i — R«; 2d ) is the Gaussian distribution for the l-th bond vector. See [5j-[7j for 
accounts of the standard Gaussian model. Can (4.29) be approximated by (4.31)? 

For N = 3,4,5, [Ajv_ 1 ] -1 / 2 is given in |86|. Based upon the latter and through 
further approximations on [Atv-i] -1 ^ 2 in (4.25) for large N, one gets, for the single open 
freely-jointed chain, results for various physical quantities: correlations among bond 
vectors, squared end-to-end distance, probability distribution for the end-to-end vector, 
behaviour of the chain under weak external stretching forces and "rubber elasticity", 
structure factor for small wave vector. Such results do agree with those implied by 
the standard Gaussian model (say, by (4.31)), which provides a check of consistency: 
see 26 for details. Further analysis on [Ajv-i] -1 ^ 2 , more detailed and considerably 
improved, have been carried out in 72-74] also for the case Mj = M (i = 1, . . . , N) and 
dj = d (j = 1, . . . , N — 1), thereby providing further confirmation of the consistency with 
the standard Gaussian model. Thus, one finds the following approximations [26, 73, 74] : 

Ri;2(JV- l)d 2 ) 



W(Ri, ...,Rjv) 

Z io) 

z 

G(N-l) 



W eq , W(R. 



N 



~k B T~ 


N-l 


2ttH 2 _ 




4-kR 3 q 


-7 


3 J 





R x ) ~ W G (R N 

d 2(N-l) 



{N/M N - 1 fl 2 
[d 3 Rjvd 3 Ri] G(N - 1 



G(R N ;Ri;N-l) 



N-l 



n ^ 



l'=2 



w, 



eq ■ 



(4.33) 

(4.34) 
(4.35) 
(4.36) 



Notice that Z and G(N — 1) are manifestly rotationally invariant. Also, the factor 
[(47ri?Q/3)] _1 cancels out a similar contribution arising from the integrations over both 



4.3 Chain with weak next-to-nearest-neighbours interaction 
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Rat and Ri in Z. 

(see 
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^R,app (T- independent) is an approximation estimate for Zr 
72 and Eqs. (C.4) and (C.2) in [73]), and it will not be relevant here. 



The representation in (4.34)- (4.36) displays quite directly the connection with the 
standard Gaussian model. One finds 26 73, 74: a) (iij • Uj) ~ for % ^ j, b) 
((Rat - Ri) 2 ) ~ (AT - l)d 2 . Eq ( |4.34| ) includes Z explicitly, even if Z = 1: such a 
redundancy is justified because (4.34)-(4.36) will help to understand the generalizations 
(8.2)- (8.3) for ds chains (|8|), in which the counterpart of Z will be 7^ 1. 



of 



Let us revisit briefly the case with external stretching forces in cluded at t he en d 
Then, Zq\ becomes, in the classical limit, Z^f , given by ( 4.23| and (4.24), 



4.1 



provided that the latter includes, in its integrand, a factor exp[— (hsT) x f J^iLi 1 ^i u i]- 



For applications, see 26 



4.3 Chain with weak next-to-nearest-neighbours interac- 
tion 

We continue to assume di = d and Mj = M, i = 1, . . . , N — 1. Here, we shall suppose 
that, besides the strong nearest-neighbour interactions through Vj = (2Bj)~ 1 Lj 2 (yj — d) 2 
with large coj in (4.1), other weak potentials (Vjj+i = v(\ yi + yi+i |)) exist between 



atoms which are next-to- nearest neighbours (nnn), v being a real function. Then, the 



total interaction potential is now: U(y) = ^jLi 1 Vj + SjLi^d Yi + Yi+i I)- The 
case in which all Kj+i be adequately strong (and, then, constraint the bond angles) 



will be analyzed in section [6j The analysis in subsections 4.1 and 4.2 can be extended 
straightforwardly to this case. We shall proceed to treat this chain in the classical limit. 
Then, (|4.23h- (14.241) hold, with Z R replaced by: 



Zr 



[dn][A 



N-l\ 



-1/2 



exp [-(k B T) 1 U nn 



(4.37) 



with U„ 



2^j = i v(d I uj + u i+ i 



Let W(r) be the end-to-end distribution 
(r = Rtv — Ri), given in (4.30), and let G(q) be its Fourier transform, q(= (0,0, q)) 
being the wavevector (compare with Hi). One has: 



G(q) 



d 3 rW(r) exp(iqr) 



/[dOKAjv-il-^exp 



N-l 
3=1 U J 



exp [-(k B T) 1 U n 



/[d^Aiv-il-^exp [-(k B T)-W n 



(4.38) 



We suppose that v is not deep but very smooth, so that its minimum times (feT) -1 is 
negative and, in absolute magnitude, smaller than unity. According to [76] , an important 
necessary condition for the replacement of quantum partition functions by classical 
ones be justified is that interparticle distances be appreciably larger than thermal 
wavelengths (" the classical limit restriction"). Let Xth = (27rh 2 /MbsT) 1 / 2 be the 
thermal wavelength. The above classical limit restriction reads Xth *C d, which holds for 
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adequate T. Let us consider the distances [(R/ — Rj) 2 ] (with I = 1, . . . , N — 2 and 
N > j > I + 1) between the / — t/i and the j — t/i atoms. For a better characterization 



of the model, the integrations in Eq. (4.37) should be performed, by assumption, over 

\2i V2 



R, 



> >Hh > °> in 



all angles consistent with the classical limit restrictions: [(R; 
general. These would seem to amount to some sort of excluded-volume effects. Then, 
(R i+2 - Ri) 2 = 2d 2 [l + Ui ■ u i+ i] > X 2 th . Then, one can reasonably conjecture that, 
for large N, Zr gets important contributions from small integration domains where: a) 
Ui ■ Uj + i is close to +1, and b) Uj • Uf+i is larger than about —1 + X 2 h /2d 2 , with, in turn, 
1 S> X 2 h /2d 2 > 0. Let pq(< 1) account for a small excluded-volume effect arising from 
the above classical limit restriction. A reasonable estimate is pg ~ 1 — X?,/(2d 2 e max ), 
with e max (0 < e max < 1) being a cutoff value for the integrations near u« • Uj+i = —1. 
Clearly, pg = 1 amounts to the absence of classical limit restrictions. We shall 
estimate A^, at room temperature (ksT ~ 0.02 electronvolts), for d ~ 5A, M ~ 12 
proton masses ( about the mass of a Carbon atom). One finds: Xth — 0.6 A. One 
has pg ~ 0.84 for e max — 0.05, while pg ~ 0.986 if e max ~ 0.5. We introduce 
p = pg exp{(kBT)~ 1 [v(2d) - v (X th /(2 1 ^ 2 d))]}(> 0). In principle, p could be larger 
or smaller than unity, but not much in both cases, due to the assumed smoothness of 
the actual v. 

The squared end-to-end distance ((xjy — x i) 2 ) f° r the chain with iV atoms can be 
obtained from G(q) upon expanding for small q through: G(q) — 1 — 6~ 1 q 2 ((xtv — 



+ 



A somewhat lengthy calculation, gives: 



((xjv - xi) 



d 1 



i+N + (l-p) N -(l + p) N 



(4.39) 



p 2p 2 (l + p) N - 2 

See appendices B and C for the derivation of Eq. ( |4.39| ) in outline. Then, for small-g: 

g 2 ((xiv - xi) 2 



G(q) ~ exp 



6 



(4.40) 



Figure 4.1 shows a comparison of G(qd)(= G(q)) as obtained from Eq. (4.40) and from 



Eq. (B.3), when we set in both p = 1. In general, as increases, there is a significant 



improvement of the agreement between formulas (B.3) and (4.40). 



Inverting Eq.(4.38), with (4.40), one finds the end-to-end distribution: 



W(r) 



27r((x 7V -x 1 ) 2 ) 



3/2 



exp 



3r 2 
2((xjv - xi) 



(4.41) 



Some interesting consistency checks of Eqs. (4.39)-(4.41 ) are: 1) for p = 1 ( by omitting 
the classical limit restrictions and assuming U nnn ~ 0), one gets ((x^r— xi) 2 ) = d 2 (N— 1) 
(random coil); 2) for p close to (that is, going to the limit of physical reliability of our 
computation and perhaps beyond, as v has been assumed to be smooth) and fixed N, 
some algebra gives ((xjv — xi) 2 ) ~ d 2 (N — l) 2 (fully stretched chain); 3) as — > +oo, 



with fixed p in < p < 1, one gets ((xjy — Xi) 



d 2 (N/p). We emphasize that 
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Figure 4.1: Comparison of G(qd), obtained from Eq. ( 4.40[ ) (dashed 
line), to that given by the exact formula in Eq. (B.3) (solid line), 
when one sets p = 1 and increasing values of N. The inset shows the 
details of the curves within the dotted region. 



when p = 1, Eq. (4.41) coincides with the well known result for the Gaussian end-to- 



end probability distribution [3j-[6] and also provides a check of consistency, a posteriori, 
of the approximations made in appendix B. On the other hand, 2) and 3) indicate 
that the limits N —> +oo and p — > cannot be interchanged. Clearly, a behavior for 
((xjv— x i) 2 ) proportional to cPN®/ 5 for arbitrarily large N (say, Flory's behavior) cannot 
be obtained analytically in the actual model with next-to-nearest neighbors interactions 
only. However, for small p and a restricted but interesting range of (large) values for 
N, Eq. (4.39) gives values which could resemble Flory's behavior numerically. In fact, 
take, for instance, p ~ 0.2 and 10 6 < N < 10 8 : then, ln[((x# - xi) 2 )/d 2 ]/(ln N) varies 
between 1.116 and 1.087. It is also interesting to compare Eq. ( 4.39| ) with the squared 
end-to-end distance ((xjy — xi) 2 )y, obtained in another model for an open linear chain, 
treated as a cooperative system |3j: the latter model in [3j is not based on Quantum 
Mechanics and is formulated directly on an analogy with the one-dimensional Ising 
model. In general, ((xjv — xi) 2 )y differs from Eq. (4.39), but in the limit N — > +oo, 
^ 2 )y ~ d 2 (N/p), that is, it agrees with the above result 3) (when, in the latter, 
= 1, v(2d) — u(A t / l /(2 1 / 2 ci)) < 0), which provides a partial check of physical 



((xjv-xij 
we set pq 



consistency for Eq. (4.39). 



4.4 Star freely-jointed polymer 



Another interesting class of macromolecules is formed by the so-called open star 
polymers [XT]. The procedures for dealing with quantum constraints discussed in 4.1 
can be generalized, without essential difficulties, and have led to a model for a three- 
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dimensional open star polymer [72], as we now summarize. We treat the case of n arms 
(n > 3), all of which start from a certain atom ("the vertex or origin of the star"), with 
position vector Ro and mass Mq. The r-th arm has N r atoms (without counting the 
"vertex"), forming a linear subchain ( 1 < r < n). The mass and the position vector 

(r) (r) 

of the i-th atom along the r-th arm are M> and R^ , 1 < i < N r ). Then, let the 
center-of-mass (CM) position vector and the relative ("bond") ones along the r-th arm 

(r) 

be denoted by Rca/> and y , respectively. One defines: 



CM 



n N r 



Mi 



lot 



M o + £E M ^' (4.42) 



tot 



--1 i=l 



and y 



(r) 



R 



(r) 



R, 



(r) 



R^i f° r (1 < * < N r ), with the understanding that, for any r: 
Ro (the position vector of "the vertex"). Below, we denote the set of atomic 



coordinates (y^) by y. The starting quantum Hamiltonian is also Eq. (3.1), with 



(r) 

obvious replacements of indices. Then, using expressions (4.42) and y , (3.1) gives 



(after various cancellations) the quantum Hamiltonian in terms of the new variables as: 
Hq,i = -(h 2 /2M tot )V* + H Q , with : 



t.2 n n N r »-2 

= -2j^<E v y <,) 2 -EE <V - v.,, ? + < -y, < i. 

u r=l r=l j=l 



with the understanding that V y ^ = 0. The total potential energy interaction 
f/(y) = V(R) is independent on Rcm- We shall assume it to be: 



f(y) = EE 



r=l i=l 



25} 



(r) 



(4.44) 



are large frequencies and eQ are the bond lengths. One employs: 



(M q 



L + 



also with M, 



o 



M Q . Also, y ? (r) =| y) 



We shall assume that the star 



polymer is in thermodynamical equilibrium at absolute temperature T. The analysis for 

(r) 

large u;- , leading to an effective description in terms of slowly- varying angular degrees of 
freedom, follows the same pattern as that in |4.1[ We assume a variational wavefunction 



N r 



having a similar structure, now with IIr=i H 
corresponding gro und-state wave function for t 



^(Vi^) tPo-{d,p), where (j)\ r> is the 

le harmonic-oscillator-like potential, 

- (r) 
similar to Eq. (4.7). 9, <p denote the set of angular variables for all y) . We perform the 



0) ■ 



variational computation and employ Peierls' inequality and so on. Finally, the classical 
Hamiltonian reads: 



2M 



E 

r=l 



Ar) (r) 



N r 



M,c 



(0^2 



) 2 + EE 



i 



=i i=i 



■IM^d. 



Ar) 
i.c 

7^) 



Ar) 



d) 



(r) 



(4.45) 
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(r) 

with the understanding that a^ r)+1 c 
star polymer is (N = Ylr=i ) : 



0. The classical partition function for the open 



k B T 
2ttH 2 



N 



ilr=ilL=i (d; 



ATM fjry 



(det£) 3 / 2 



dQ] [det(i4 r) (B 



-1/2 



(4.46) 



The elements of the actual (N) x (N) matrix B are readily obtained from Eq. 

(4.45), and we shall omit them. The internal energy U of the star polymer is obtained 
through procedures similar to those in the previous section. One again finds energy 
equipartition: U = NksT. 
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Chapter 5 



Single-stranded closed-ring 
freely-jointed chain: getting 
constant quantum zero-point 
energies of hard degrees of 
freedom 



The relative simplicity of open freely-jointed chains, summarized in 4.1, may cause 
the impression that the generalizations of those methods to interesting quantum 
macromolecules, with other kinds of constraints, be always easy or direct. Unfortunately, 
such an impression is not correct. We shall now proceed to another kind of constraint in 
quantum macromolecules, which poses new and considerable difficulties and which will 
have to be tackled, in order that those generalizations could be accomplished, at least 
up to certain degree. 

We shall construct a quantum-mechanical model for a three-dimensional closed- 
ring (freely-jointed) molecular chain. To fix the ideas, we shall suppose that the 




Figure 5.1: Single-stranded closed-ring macromolecular chain. 
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atoms in the chain and/or the bond lenghts are not all identical to one another. Le t 

1/2 



VN 



denote the distance between the atoms 1 and N (see figure 5.11). 
Notice that j/jv is not an independent variable. A first warning of the difficulties to be 
met is that all variables yi {I = 1, . . . ,N — 1), which are independent, and yiy, which 
is not, should be treated on the same footing. That is, even for non-identical atoms 
and/or bond lengths, the choice of where the chain was closed to become a ring should 
be completely arbitrary and physically irrelevant. We shall choose the Hamiltonian 



Hq to be also given in Eq. (3.2). If our purpose is to extend the quantum-mechanical 
variational computation of section 4 so that we end up with a quantum model for a 
closed-ring chain: what should U(y) be? We shall restrict to harmonic-oscillator-like 



potentials. It seems reasonable that U(y) now includes the right-hand-side of Eq. (4.1) 
plus an interaction potential between the atoms 1 and N, namely, {2B^)^ l uj 2 N {yN — dN) 2 
( B N = Aftf 1 + Mf 1 ). d N and ojn are another bond length and a large frequency, 
respectively. Another question is: What should the trial variational function <p{y) (the 
counterpart of ( |4.4| )) be ? A reasonable variational wave function <fi{y) for a chain, 
which should become a closed ring in the oji — > oo limit (i = 1, .., N — 1, N), could be 
chosen as: 



N 

0(2/1,2/2, ■ ■ - ,2/iv) = n^( y; ) • 
i=i 



(5.1) 



Thus, both U(y) and (j>(yi, 2/2, • • • , Un) include all variables yi (I = 1,...,N — 1,N) 
on the same footing. Each (fti(yi) is chosen as the ground state wave function for the 



corresponding harmonic-oscillator-like potential see Eq. (4.7). Then, the (normalized) 
total trial radial-angular wave function of the chain reads: 



4 = *(y) 



iV 



.1=1 



(5.2) 



The normalization condition f[dy] \ $(y) | 2 = 1 ( with [dy] = Y\iL\ l d^yi) becomes in 
the limit uj\ — > 00 (by virtue of Eq. fO 



dn]|w(^)| : 



8(vn - d N ) _ 



d 2 

a N 



(5.3) 



[dO] being given in (2.13)). Except for (5.3)), ipu{0, ip) is fully arbitrary. The restriction 
given by 5(yiy—diy) (after having imposed those related to I = 1, . . . , N— 1) is interpreted 
as: 



VN = VN(0,ip) 



AT-l 

/I didjUiUj 



1/2 



d N 



(5.4) 



which establishes the closed-ring constraint among all 2(N — 1) angles contained in 8, (p. 
This displays, at this level, the difficulties of the quantum closed-ring chain. 
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Our first aim is to evaluate, as all frequencies i = 1, .., N — 1, N grow very 



large, (fc.iZgS) = f[dy]$(y)*H Q $(y), by using Eqs. (3^, (48), (5^ and (5^. 
After lenthy computations and cancellations, one finds (see appendix D in this review 
and (71]): 

N-l 



i=l 



(a* • ai) 



S(vn ~ d N ) | 



J2 

U 7V 



iV-1 



1 



[dn]<(^) 



8=2 



didi-i 



$(yN - d N ) 



d 2 

a N 



(5.5) 



for any ij) a {6,(p). a; are the same as in Eq. (3.4). Here, {0^(h)j denotes the 

total contribution of the set of all terms which do not depend on the frequencies 
and which do include, at least, one operator acting on the delta function: they and 
some additional calculational details are given in appendix A in [71]. Therefore, all 

contributions contained in ( 0^ c \h)^j are proportional to, at least, the first power of H. 

(not employed for similar contributions 0(H) 

The meaning of, say, 



We use here the specific notation 



for other chains) in order to facilitate comparisons with 71 
[(aj • aj)Vv(0) <p)] fi(UN — d^ in Eq. (5.5) is the following. First, the differential operator 
(a^ ■ a^) acts upon ip a by regarding all angular variables 0\, <pi, I = 1 , . ., N — 1, as if 
they were independent on one another (that is, as if the constraint 5(un — disf) were 
not operative). After the operator has acted upon ip a with that understanding, then 
8(vn — d^ acts and implies that there is one relationship among those 2(N — 1) angles. 
A similar interpretation applies for [(aj_i ■ ai)i() a (0, cp)] 5(dn — d/v)- These operations 
are genuine consequences of the closed-ring constraint at the quantum level. 

It is crucial to notice that all dependences of (<£, Hq&) on the frequencies appear only 
in the constant Eq = Yli=x ^rS m the right hand side of Eq. (5.5). It is rewarding that, 
after the computation and the cancellations (outlined in appendix D), Eq is independent 
on angles, equals the sum of all zero-point energies associated to all harmonic-oscillator- 
like potentials (including the N-th. ), and that Eq depends on all frequencies Wj, 
i = 1, .,..N — 1, N on the same footing. The relevance of all that is explained in the 
following paragraph. 

A priori, one could have expected that, due to 5(yN~dN), the calculation of (<J>, Hq&) 
would give rise to terms linear in the frequencies multiplied by integrals containing 
functions D{ depending on, at least, some of the angles 6i, . . . , 9n-i, <Pi, • ■ ■ , fN-i ( sa Y) 
to pieces like 

oji f[dn}4>*(6,(p)Di(8,p)ip a (e,p), with % = l,..,N - 1,N). Then, these integral 
terms would have implied wild variations of the integrand in J[dfi], as they would 
become amplified by the diverging vibrational frequencies. This would upset the 
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reliability of the result for (&,Hq<&), as commented in 54, 56: see also 26 
and references therein. However, if, by virtue of the algebra and cancellations in 
the variational computation, it would turn out that all functions D t be constant 
(independent on the angles), then, the result of our calculation for (<&,Hq&) would be 
reliable, and we could separate, in a physically meaningful way, the very large vibrational 
frequency contributions (which would become constant, by virtue of the normalization 



condition in Eq. ( 5.3 ) ) from the smaller rotational contributions (which do imply angular 



variations). The important result displayed by Eq. ( 5.5 ) is that, in fact, all terms linear in 
the frequencies go multiplied with factors which are, indeed, constant, as a consequence 
of our choice of the symmetric interaction U(y), the variational trial wave function in 
(5.2), the algebra, the cancellations and Eq. ( |5.3[ ) [71] . This solves the difficulties a), b) 
and c) at the end of subsection 2.2. 



Let us make use of the variational inequality (3.7). Then, the quantum partition 
function Zq for the actual closed-ring chain fulfills: 



Z Q > exp [-(^r)- 1 ^] ■ Z { 



(c) 



z. 



Q 



E>xp 



k B T 



(5.6) 
(5.7) 



The quantum angular Hamiltonian Hq for the closed-ring chain reads: 



H. 



Q 



N-l 

E 

i=l 



Bi 
2 



(gj • aj) 
d 2 



N-l 

E 

i=2 



Mi 



did;, 



(5.8) 



with the operational meaning for (a^ • a,) and (a^_i • a,) indicated above. We emphasize 

(c) 

that Zq can be regarded as the three-dimensional angular quantum partition function 
for the closed-ring chain. A posteriori and by grouping the resulting terms in Eq. (5.5), 



the latter can be reformulated. After some lengthy algebra, Eq. (5.5) becomes 751: 



($,# Q $) = E + O^(h)x 



f 1 at— l „ 

i,j=l ■> 



5{un - djv) 



d 2 

a N 



(5.9) 



(with 0^{h) i — > as h — > 0). From (|5.9[), invariance under overall rotations in three- 



dimensional space can be shown to hold for Zq\ See 

We have also carried out (as an unpublished work) the extension for a quantum 
closed-ring freely-jointed chain, by using Morse potentials, instead of harmonic- 
oscillator-like ones. Then, U(y) now equals the right-hand-side of Eq. (4.3) plus 
another Morse potential between the atoms 1 and N, namely, Vm{un)- The resulting 
computations and cancellations follow the same pattern as those for harmonic-oscillator- 



75 



like potentials (although they are more complicated). Eqs. (5.5)- (5.8) continue to hold 



now with the constant Eq standing for YliLi ^M,i 



n=0- 
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M in the chain, i = 1, . . . , N) and, 
d, j = 1 , . . . , N) . The comments 
on indistinguishability in 4.1 keep their validity here. A new deep issue related to 



Let all atoms be identical to one another (Mj 
for simplicity, let all bond lengths be equal (dj 



indistinguishability appears here, arising from the fact that any atom can be considered 
as the first one. In order to deal with it, an overall prefactor (scaling as 1/N, eventually) 
could be introduced in the quantum partition function, which would not alter the 
variational computation. 



Eqs. (5.7) and (5.8) solve the difficulties a), b) and c) at the end of 2.2 for the actual 



chain (which justifies their derivation), but their complicated structures make their 



practical use very difficult. Eq. (5.7) will become considerably reduced, if one proceeds 



to the classical limit approximation for the angular degrees of freedom (by assuming 
Eq. (4.19) for I = 1,..,N). Notice that (O^ \h)) vanishes in that limit. The quantum 



■(c) 

Hamiltonian Hq for the closed-ring chain can be approximated, in the classical limit, 
by the same H, 



(o) 



r(SJ 



in Eq. (4.20) [with the same classical aj 

tiplied by 5(y N 



in Eq. (4.21)] as for the 



dN)/d 2 N . An analysis of 



open chain: of course, He ' goes now mu 
this classical limit has been carried out in 

be approximated by the following three-dimensional classical partition function for the 
closed-ring chain: 



7lj , and we shall quote the result, zj^ can 



J C 



2irh 2 



N-l 



N-l 
1 



(det5)3/2 



dfi 



S(vn - d N ) 
d 2 

a N 



[A 



N-l\ 



-1/2 



(5.10) 



with the same (det£) 3 / 2 and [A 

N-l] as in Eqs. d4.26b and d4.27|. 5(y 



-In - d N )/d 2 N 

appeared in the exponent in the quantum Eq. ( |5.7[ ), but it is located downstairs in the 
classical approximation. 

We shall comment briefly about an important necessary condition for approximating 
quantum partition functions by classical ones [76]: interparticle distances should be 
app reciably larger than thermal wavelengths (" the classical limit restriction", like in 
4.3 ). Let us employ the distances [(R; — R-j) 2 ] ^ 2 between the l—th and the j—th atoms, 
which are to be expressed in terms of ( Eq. (1.2) ). In general, the integrations in Eq. 



(5.10) should be performed over all angles consistent with the classical limit restriction: 

1 /2 

all the atomic distances [(R; — Rj) 2 ] should be appreciably larger than all thermal 
wavelengths \h,j = (2irh 2 /MjkBT) 1 / 2 , j = 1, .., N. The approximation of Z^q by Z^q is 
justified only when they hold. An estimate of thermal wavelengths at room temperature 
is given 4.3 Then, those restrictions would imply the exclusions of some (eventually 



rather small) angular domains in the integrations in ( 5.10| ). See 71 and [72| . 

Eq. (5.10) yields energy equipartition: U = (N — l)ksT, as for the open chain. 



Topological constraints and knots play an important role in the properties of 
closed-ring chains 187, 88 . We would expect that the variational quantum- mechanical 



(c) 

computation keeps its validity regardless topological constraints, and that Zq would 
count and include all kinds of topological conformations (knots) that a closed-ring chain 
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may adopt in space, and so forth for Zq (under the conditions of validity of the classical 
approximation). Such studies lie outside our scope, as commented in [71~]. 



Chapter 6 

Single-stranded open 
freely-rotating chain: getting 
constant quantum zero-point 
energies of hard degrees of 
freedom 



The Hamiltonian Hq continues to be given in Eq. (3.2). We shall suppose that, not 
only nearest-neighbour atoms interact through stiff harmonic-oscillator-like potentials 
Vj = (2Bj) oj 2 (yj — dj) 2 {yj =\ yj |) with large vibrational frequencies C0j and lengths 
dj(> 0) [75], but also that similar potentials exist between atoms which are next-to- 
nearest neighbours: Vjj+i = 2 _1 Bjj + \uj 2 Yj + Yj+i \ ~dj,j+i) 2 - Here, are 
other frequencies, djj + ± are lengths such that | dj — dj + \ |< djj+i < dj + dj + \ and 
Bjj + i(M^ 1 + M~ +l ) = 1. Then, we suppose: 

N-l N-2 

For suitably large Ujj+i, Xl^^i 2 ^j'J+i hinders part of the allowed internal rotations in 
the macromolecular chain, which becomes a freely-rotating one. On physical grounds, 
^2!j=i could approximate for the effect of the covalent bonding due to successive 

single pairs of shared electrons, which produce precisely those hindrances [3}|5]. Weak 
next-to-nearest neighbour interactions were treated in |4.3| , while strong ones are 
considered in this section. 
We shall suppose: 

fiw, » k B T , hui > , (6.2) 
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Constrained 
y : :; ' length 

#yj+yj+i 



Figure 6.1: Freely-rotating single-stranded macromolecular chain. 

h 2 B l 



hwi i+i > k B T , fojju+x > 



4 



(6.3) 



for any I. Physically, one would expect that hw\ be larger than huiu + i, because the 
interaction responsible for the existence of bond lengths should be certainly stronger 



than that giving rise to constrained bond angles. Based upon (3.6), (6.1) (6.2) and 



(6.3), the following variational computation will lead to a model for an open freely- 
rotating chain, in which all dj and djj+x be given constants. We choose the variational 



(radial-angular) wave function <&(y) as: 

$(y) = (pnniy^nnniy^aiO,^) , finniv) = J| 



N-l 



(6.4) 



1=1 



4>i(yi) being given in (4.7). In turn, we choose: 



y) = rL=i 2 0M+i- Each 01,1+1 is 



also chosen as the real (Gaussian) wave function associated to Vij+x'- 

^1,1+iB^i+i 



1,1+1 



^1,1+1^1,1+1 


1/4 


exp 





2h 



yi + yi+i I -d, 



1,1+1, 



(6.5) 



As all frequencies uji and w^+i become suitably large, Gaussians approach Dirac delta 
functions (w rad = [JJ^ 1 <^ 2 <%z - d t )]): 



w. 



ang 



I H J nn (y)4>nnn(y)\ ~+ 
JV-2 

Y[ S (\diui + dj| + iU| + i| - dij+i) . 



(6.6) 
(6.7) 



l=i 



5(| diui + di + \U[ + i | — dij + i) is equivalent to constraining the bond angle between u/ and 



u/_|_i: see Eq. (6.15). The complex functions t/j a (6,(p) (a now being a set of quantum 
numbers) are arbitrary, except that: i) they are periodic in each <pi with period 2tt and 
independent on any yi, and ii) they are normalized through the scalar product: 



[dfl]w ang ipi(6,ip)*ip 2 (9,ip) . 



(6i 
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One evaluates, as all frequencies uji i = 1, . . . , N — 1, and / = 1, ••, N — 2, grow very 

large, the quantum expectation value (§,Hq<&) = f[dy]$(y)*HQ$(y), by using Eqs. 
(6.4)-(6.8). The computations and cancellations [75], which are quite lengthy, follow 



the same pattern as those in |71| and embody similar cancellations, now applying the 



angular constraints in (6.7). One novel aspect is that, by the end of the computations, 
one manages to group the resulting terms into some, quite compact, structures. The 



result is 75 



($, H Q <5>) =E + (W, H% r) ip„) + C)W(fi) , 



(6.9) 
6.10) 



for any normalized if) a (6,(p) fulfilling the above requirements. The constant Eq equals 
1~ X ^i + ^La=\ 2~ 1 frcoi ! i+i (namely, the sum of the zero-point energies associated 
to all Vi and Vjj+i). In turn, the angular Hamiltonian Hq r ' is given through: 



w, 



angHq 



(fr) 



N-l 



It, 



Bji_ 



(Wangej)} 



3.11) 



with = ihui — aj, I = 1, . . . , N — 1. 0^ r '{K) denotes the total contribution of the set 
of all remaining terms which do not depend on any of the frequencies u>qi and Woii+l- 
0^ T \h) is proportional to ft? and it does not contain differential operators acting upon 
ip a (9,(p) (that is, w ang O ang {ft) acts multiplicatively on ip a (9,(p)). The explicit form of 
0^ T \ft) , which is not relevant here, is given in appendix A in 



75 



The structures in 



(6.9)- ( |6.11 ) solve the difficulties a), b) and c) at the end of subsection 2.2. and ressemble 



formally those which appear formally for the freely-jointed chain (recall (4.13)- (4.15) 
), which appears to be rewarding. Notice also the formal analogy with (5.9). We 



emphasize the crucial importance of having chosen the variational trial wave functions 



adequately (through (6.4), (4.4), <j) nnn (y) and (6.5)). Otherwise, one may obtain either 
contributions depending on the unconstrained soft variables or results which, even if 
independent on the latter variables, do not yield the physically expected quantum zero- 
point energies. An example of the latter situation was given in [26] (subsection 6.2)), 
where a variational trial wave function containing factors differing from ( |6.5[ )) was used. 
We continue to apply Eq. (3.7). In our case, $o- = $(y) and (<£ CT , Hq<&o) are given 



in Eqs. (6.4) and (6.9). 



Then: 



Zq > e W [-(k B T)- L E ] ■ Z Q J 



4> /r> - E 



exp 



(kBT)-^ / ldfl]rA0^)WangH { o fr) M0^)+O {fr) (ft)} 



3.12) 



(Q.13) 



where Zq can be regarded as the effective three-dimensional quantum partition 
function for the slow motions of the unconstrained angular degrees of freedom of the 
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freezing chain. <> is a Haitian operator (with respect to the scalar product 



inEq. (6.8)): 



ifr) 



(Vi,^ /r V 2 ) 



V>2) = (H, 
N-l 



ifr) 



ipl,ip2), for any ipi,ip2- One has: 



(6.14) 

a hi , 

i,3=l 

It can be justified that and, hence, (with O^(h) discarded) are invariant 

under overall rotations in three-dimensional space. The above properties have indeed 
been established, at the price of performing various partial integrations and handling 
carefully w an g (as the latter involves Dirac delta functions). See appendix B in 75 
We shall suppose that all angular wave functions ij) a 



, tp) are the complete set of a 
orthonormalized eigenfunctions of Hq^ so that, upon discarding {fr \h), Eq. (4.18k 

'It[e3cp[-(JbflT)- 1 fl'^ r) 



becomes Z 



ifr) 



I) 



jq — j-l [c^jj L — \rv B ^ j J.J.Q j. It is convenient to express the 2(N 
angular variables in and <p in terms of another set of 2(iV — 1) more suitable 
The latter will be chosen to be: 9\,... , 0jV-i> fo{= (N — 1) _1 Ylj=i fj) 



ones. 



and 



UoU 



N- 

r- 



cos Oj cos 9j + i + sin Oj sin Oj cos(ipj + i — <fj)), j = 1, .., N — 2. One has: 



, jj+i]d^o[n^=i 1 the Jacobian J depending on all Oj and Pjj+i- 



0- 

The angular constraint reads ((u.jUj-+i)(°) = (2djdj+i)~ 1 ((^ 3 - +1 — 



iV-2 

n 



dj,j+i 
djd j+1 



N-2 

n 

3=1 



(Uj 



u, + i) (0) ) 



3.15) 



Specifically, the constrained bond angles are all those given by cos -1 [(u.j • Uj+i)(°J], 
j = 1,...,N — 2. Under (6.2) and (6.3), we shall study the transition to the 
classical limit, which is simpler if one has got rid of all f3jj + i, previously Then, we 



start from Z, 



ifr) 

Q 



Tr[exp[— {ksT)^ 1 Hq T ^] : apply to it Peierls' variational inequality 
(3.7), employing now the complete set of all ip a (0, <p)'s depending on £q, . . . , 0jv-i and 
ifo (but not on /3j,j+i) and use Eq. (6.15), for integrating over and getting rid of 
all Pjj+i- The resulting quantum partition function has an adequate form for taking 
h — > 0. About room temperature, the individual slowly-varying internal rotations about 
bonds, which remain unconstrained, have typical energies which can be estimated to be 
smaller than the vibrational energies huii and fuoij+i (3JJ6, 68 . Then, we shall suppose 
that ksT 3> (h 2 Bi)/(df), so that an appreciable number of excited states for those 
unconstrained rotations may be occupied and quantum operators and statistics can be 
approximated by classical ones. Accordingly, all quantities of order h or higher (like 
0^ T \K) ) can be neglected directly. We shall omit the analysis, which is given in [75] . 
One finds the following classical (c) Hamiltonian: 

flW = I Y^aW-aW. (6.16) 



h3~- 



" n ij ifr) 



3,c 



a'^ s are classical variables ( arising from the classical limit of ej): a, 
{{N — 1) sin#j)~ UjjjP^]. Pg i and P vo are the classical momenta canonically conjugate 



ifr) 
i,c 



>8M. + 
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to Oi, i = 1, . . . , N — 1, and ipo- in the classical limit, one gets the classical partition 
function Zr^ r \ 



Z (fr) 



(2Tth) N 



N-1 



d9 i dp e, 



1=1 



difodPpQ exp 



H, 



3.17) 



One performs all Gaussian integrations in Z, 
and over ipo . Then: 



(fr) 



over the classical momenta Pg i and P, 



Z (fr) 



1 



(2irh) 



N 



2ir[2Trk B T] 



N/2 



N-1 

n d0 i 

. i=i 



[D 



N 



-1/2 



3.18) 



Djv is the N x N symmetric matrix formed by the coefficients of P Vo and Pf 



1, ...,JV- 1 in il. 



(fr) 



The internal energy C7 of the classical freely-rotating chain 



can be obtained directly from (6.18) ( like for the classical freely-jointed chain). One 



finds readily energy equipartition, namely: U = (N/2)KbT 26 72 



A simple and interesting example of an open freely-rotating chain is polyethylene 
(the synthetic polymer ... — CH2 — CH2 — CH2 — CH2 — ...., each CH2 being regarded, for 
simplicity, as one of the "atoms" ) . We accept that bond lengths do not vary appreciably 
along the chain [dj ~ d) and that the same holds for bond angles. For polyethylene, the 
bond length is d = 1.54 x 10 _1 nanometers and the bond angle cos -1 Pj j+i is about 70 
degrees and 32 minutes (for any j) |6j. 

Let dj = d, j = 1, ..,N — 1, and Mi = Mo, i = 1, ..,N. It now follows (using some 



long-distance approximations ) that Eq. (6.18) implies the existence of a persistence 
length dpi (> d): see |26| (subsection 6.4) for a qualitative estimate. So, on a suitably 
large length scale, the classical freely-rotating chain described by Eq. (6.16) can be 



approximated by a freely-jointed one having N p i bonds (N p i < N — 1) and a new bond 
length d p i. 



The procedures in this section leading to Eqs. (6.9), (6.10) and (6.11) have been 



generalized, in outline, to single-stranded closed-ring freely-rotating chains and to single- 
stranded open freely-rotating chains with further constraints (helical-like, star-like...), 
all at the quantum level [75] . We remark that one always finds structures similar to those 
in Eqs. (6.9), (6.10) and (6.11), with different w ang associated to the corresponding 



constraints other than those for constrained bond lengths. That is, the fact that bond 
vectors are constrained is not represented by any contribution to w ang , but by the very 



structure 2 1 Yld,j=i{ B ij / d i d j) e i 



jointed case, if w ang = 1: compare with. (4.15). Eq. (5.9) was one particular case 
of (6.10). As an example, we remind that the DNA of the bacteriophage (j) X 174 is 



ej in (|6.11|). Thus, (|6.11|) also includes the freely- 



single-stranded and forms a ring iTU (and, of course, it has angle constraints, so that to 
regard it as freely-jointed chain would be too crude). 



6. Single-stranded open freely-rotating chain: getting constant quantum 
44 zero-point energies of hard degrees of freedom 



Chapter 7 

Double-stranded open chain 



Double-stranded (ds) macromolecular chains (specifically, dsDNA) have an absolutely 
crucial importance in Molecular Biology [2j [3j [5j [lOj [l2j. This fact has given rise 
to an enormous and permanent research activity on Chemical Physics and Physics 
of dsDNA: see, for instance |15l 28, 89-93 and references therein. Motivated by 



dsDNA, we shall undertake in this section, in outline, a study of double-stranded open 
macromolecular chains, by extending the methods of the previous sections. We shall 
apply its consequences to open dsDNA, in [8} 



7.1 Double-stranded open quantum chain: Some general 
aspects 

We shall consider a model for a single ds open macromolecular chain, also based upon 
Quantum Mechanics. Each of the two strands or individual chains has N atoms. The 
mass and the position vector of the z-th atom (1 < i < N) in the r-th strand (r = 1, 2) 

(r) (r) 

are M- and R- . We start from the quantum Hamiltonian operator: 



He 



r, 2 



N 



EE 

r=l i=l 



M 



(r) 



V 



2 



u, 



r=l 



(r) 



r=l 



(r) 

Ut is the total potential energy among neighbouring atoms in the r — th strand. It can 



be chosen as either Eq. (4.1) ( harmonic-oscillator-like potentials) or Eq. (4.3) ( Morse 
potentials). We recall that the Morse potential gives a qualitatively adequate effective 

(r) 

interaction between two nucleotides. Us. is the total potential ( due to next-to- nearest 
neighbours interactions) accounting for the most important angular constraints in the 
r — th strand: 

N-2 



/ (r) i 



(7.2) 



shall assume that: i) \ which is the strongest potential (with the largest absolute 



This Kf,j+i can be either weak, like in 4.3 or strong (like that in Eq. (6.1)). We 
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magnitude) in (7.1 )), will give rise to constant bond lengths (yj 

W /i (r) (r) (r) (r) 

neighbour potentials V-j +1 (\dj Uj ' •' - 
strong, with a unique deep minimum 



n) the nearest- 



+ c?yYi u }+il) can be either weak or moderately 



(r 0) (r) 

value p) 1 '"'. V?> is a residual interaction in the r-th strand, not included in either of 



so as to constrain uV 



the stronger interactions C/^ or Ua 1 (and weaker than Ua 



(r) 



= F ds (Rf ) ,..,R^,R^,..,R^) is the total potential energy among atoms 
belonging to different strands: it is supposed to depend on the relative distances 
thereof, so that it displays overall rotational invariance. On physical grounds, Vds should 
approach zero as all relative distances between any atom of one strand and any other 
atom in the other strand become very large. This condition is fulfilled if one takes Vd s 
to be a sum of Morse potentials. Physically, the physically relevant values of Vd s are 
certainly smaller than those in Ylr=i + Ylr=i Ua ■ 

Let the overall center-of-mass (CM) position vector of the ds chain and the center-of- 
mass (CM) of the r-th strand be denoted by Rcm and R^M) respectively. The relative 



.(1) „(2) 



,(2) 



(r) 

position vectors between both R CM is y 



r( 2 ) _rc 

^CM ^CM 

(r) 

vectors along the r-th strand are denoted by y,- , j = 1, ..,N 



,(i) 



The relative ( "bond" ) position 
1. One has: 



td Z^r=l hn=l M i K i M 
ti-CM = TT , M tot 



Mtot 



R 



(r) _E^mTr 



CM 



2 N 



r=l i=l 



R i+1 - ilj 



(7.3) 
(7.4) 



(r) 

formed by all R- , 1 < i < N, r = 1,2, we shall find also interesting the new set of 

(r) 

independent position vectors formed by Rcm, all yj , j = 1, .., N — 1, r = 1, 2 and y. 
Then, by using (7.3) and (7.4), Eq. (7.1) can be recast, in terms of the variables of the 



to be compared to Eqs. (1.1) and (1.2) for a ss chain. Having started with the set 




Figure 7.1: Double-stranded open macromolecular chain. 



7.2 Radial variational computation 
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\2 



new set, as: H QA = -(^/2M toJ )V^ CM + Hq, with: 

*« ^ -^(v y ) 2 -E^(v y ,,) 2 -EE^(v y ,.-v ya , 

r=l r=l i=2 

- E ^(v y « / + E c/ 6 (r) + E + E y l r) + ^ . ( 7 - 5 ) 

r=l r=l r=l r=l 

( M red = Y J 2 r=i{T.i=i M l T) y 1 )- Here, the contribution -2i|^(V y ) 2 , represents the 
kinetic energy associated to the relative motion of the centers of mass of both chains. 

We shall suppose that the ds chain is in thermal equilibrium, at absolute temperature 
T in an interval about and not far from 300 K), and that l/jf^ = ^iLi 1 ^Af(j/i )> Vm 

(r) (r) 

being the Morse potential with parameters D - , a - and uj'-. We shall consider typical 

orders of magnitude, which may apply for dsDNA, at least qualitatively. We suppose 
(r) 

that all Dy are of similar orders of magnitude and that the same holds for all uj r - 

and for all v^ +l . max (= Max| V^^J). We assume that Dy are appreciably larger than 

ah v jj+i- ma x ana - that, in turn, all Vj r j +lmax are larger than all tujj . Physically, 
the interaction responsible for the existence of constant bond lengths should always 
be adequately stronger than that giving rise to constrained bond angles (recall the 
comments in|6]). Moreover, we focus on the case where &bT is less or, at most, of the 
order of hujj V ' and 

(r) , , 

tkjf' > . 7.6 

3 (4 r) ) 2 

Typical values employed in various analysis of macromolecules and of DNA are consistent 
with the above assumptions for T ~ 300 K and somewhat above (3j [6j [9j [93]. For 
instance (with 0.6 kcal/mol ~ 0.025 eV, corresponding to T ~ 300 K): Dj about 100 

kcal/mol (or larger), cfp ~ 1 to a few A, orp ~ a few A -1 , htd^p about 3-10 kcal/mol (1 

On the other 



93 



nanometer = 10A) and atomic masses typical in Organic Chemistry 

(r) 

hand, all Vj - +l . max are supposed to be somewhere between 12 and some value smaller 
than about 100 kcal/mol (6J. 



7.2 Radial variational computation 



We have performed two radial variational computations. The (normalized) trial 
variational total (radial-angular) wave function of the open ds chain reads: 

2 N-l 



n n 



r=l 1=1 



Vv(y; • 



(7.7) 



?• , tpy are the angles of y- . 9,ip denote the set of a 



variational computation 
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ft , ipf~^ . In the first radial 



are similar to Eq. (4.7) (harmonic-oscillator- like). In 
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the second radial variational computation |73|, <ft^ correspond to Morse potentials. The 
normalization condition f[dy] | <3? | 2 = 1 (with [dy] 
in the limit in which ( |4.10[ ) or (4.8) hold: 



dhiUUUl'i'd^) becomes 



hi 



dSl]\My;9,<p)\ : 



(7.8) 



where [dfi] = Yit=i llz^i 1 d(p^ doi smO^' 1 . The first radial computation of (Q,HqQ) 

(r) (r) 

(when all vibrational frequencies u\ of the harmonic-oscillator-like potentials Ui 



,(r) 



grow) is similar to that leading to Eq. (4.13). The second radial computation 173 



of ($, Hq<&), using Morse potentials, has been more difficult to accomplish. In both 



radial variational computations, one applies (3.7) for the quantum partition function 



of the ds macromolecular chain, Zq, determined by Eq. (7.5). All these lead to 72 73 : 



Z Q > expH^T)- 1 ^] • Z, 



7 (ds) 



(ds) 



£exp J [dO] J d 3 yrAy-,0,P)H { Q S) My-,9,<p) + OW(h) 



(7.9) 
(7.10) 



Here, Eq denotes Y^r=i^A=i^ ~^T~ ^ or ^he first (harmonic-oscillator-like) radial 
variational computation. For the second radial variational computation, we have that 
^0 = Y?r=l Yh=i E M,i, n =o> witn E M,i,n=o standing for the energy (< 0) of the lowest 
(n = 0) bound state for the Morse potential, constraining the i-th bond length in the 
r-th strand. The quantum angular Hamiltonian Hq S \ which is exactly the same for 
both radial variational computations, reads: 



(r) 



H, 



k2 2 pH 2 N ~ 2 dW 
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(r) 



r=l ^d 1 
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r=l i=2 



(r) 



r=l 



2d 



(r) 
N-l 



r=l 



r=l 



,(r) 

'j, 

j(r) 



a (r) • 



(0 



(7.11) 



The operators e- are analogous to those in Eq. d4.15b. X)r=i ^ 



Eq. (|7.11[) is the restriction of ^ + Er=i ^ + V ds in Eq. (|7.5[), when y\ T> = d 



r(r) 



+ E'=i^i (r) + ^s m 



for any r = 1, 2 and I = 1, . . . , N — 1, the cQ s being the bond lenghts. The remainder 
0( ds \(h) is the set of all remaining contributions, which do not depend on the frequencies 
and are proportional to some positive power of K: 0^ ds \{K) — > as h — > 0, and we shall 
disregard it in what follows. By taking the ^ CT (y; 0, y?)'s as the complete set of all 



i(r) 



orthonormal eigenfunctions of Hq' , one finds: 



,(<fe) 



Tr[exp[-(k B T)- L H. 



r( ds )l 



(7.12) 



Sec 73 for the case in which more bound states are employed in the second radial 



variational computation. Zq 8 ^ can be regarded as the effective quantum partition 
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function for the three-dimensional double-stranded chain, in terms of slowly-varying 
degrees of freedom, provided that those associated to bond angles could be treated as 



soft variables (say, like in 4.3 ). We shall be concerned with bond angles as either soft 
or hard variables in 17.21 



(r) 

Let four external stretching forces f{ and f 



and Rj^, respectively. We suppose that Ylt= 

x (ff ; + f# ; ) = 0, so that there is no net 



N ' 



(r) , ,(r)> 



1,2, act upon the atoms at R^ 



(r) 



force upon the overall CM. Then, Hq in (7.5) is replaced by: 

2 JV-1 



Hr 



F (2) 



f(2)> 
£ JV 



(7.13) 



=1 i=l 



Here, we have employed: R^ = R(r) 



I. " " /. 



M^ r) and 



The quantum-mechanical variational computation goes through as above and yields a 
new quantum partition function: Zgf = Tr[exp[— (ksT^ 1 Hq S ^]], with: 



(r) 



, ^N-l „,(V)„(r) 
'CM 2^i=? 

1, .., i — 1 and j = i, .., iV 



a^/yj where (£ 



JV 



M, 



s=l 

1, respectively 
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2 JV-1 



W„(r) 
l,i 



W u« + (ff)+i 2 ))y. 



=i i=i 



(7.14) 



7.3 Constraining bond angles and classical limit 

First, under the assumption that X^r=i be weak like in 4.3 ( the bond angles being 
regarded as soft variables ), let us turn to the classical limit. Then, w e tre at the variables 



(r) 

e\ , 9, (p in each strand like in the freely-jointed chain (subsection 4.2). One can also 
proceed to the classics, 
object, for large N 



72 



limit for y and —iKS7 y , as each individual chain is a very massive 



Then, Hq S ^ in (7.11) becomes, in the classical limit: 



H (ds) 
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W\2 \ x ^ ^ 
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r=l i=2 za j 
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(r) 



a i,c 



_(r) 
d i+l,c 



(r) 
JV 



(7.15) 



(r) 

The a- s like in Eq. (4.21), with components P( r ),P 



t.c 



corresponds to the classical 



<p. 



(r). 



The classical vector 717 



limit of —iKVv- We shall omit the associated classical 



partition function, which can b e wr itten by extending straightforwardly ( it c an be 



recovered from Eqs. (|7.16|) and (|7.17|) below, provided that the function Fj 



(■<•> 



be proportional to exp[-(£; B T)- V^Gd^ + 4+l u 5+ll 



in (7.17) 



Second, we now come to the physically very important case in which E^=i f/i (even 
if weaker than Ylr=i ■> constraining bond vectors) are still adequately strong so as 
to constrain the bond angles (hard variables, as well). Thus, let Ylr=i Ua"^ be a sum of 
harmonic-oscillator-like potentials. The radial- angular quantum- mechanical variational 
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computation for double-stranded open freely-rotating chains, in which the degrees of 



freedom for J2 r =i ^6 



(r) 



+ X)r=l ^» are treated quantum-mechanically both on the same 



footing, has been carried out in 75 , in outline. Alternatively, Y%=i Ua"^ could be a sum 



of other Morse potentials having sharp minimum for |yj- 



Mi 



dfanduf .ufjj 



P. 



WO. 



the corresponding radial-angular variational computation has been outlined in [73], with 
results similar to those in |75| . 

The last step in the second radial- angular variational computations with bond angles 
as hard variables is to obtain, upon performing the transition to the classical limit, a 
classical partition function which could be handled without unsurmountable difficulties. 



This is more difficult than in the transition yielding Eq. (7.15). We shall limit ourselves 

for 



to the case /3^ r '°' > close to +1 for any r and j, to be assumed below: /3j r '°^ 



0. 



any r and j corresponds approximately to B-DNA |3|). Then, we arrive, as a result of 
the radial-angular variational computation, at a reasonable classical partition function, 
given below in Eqs. (7.16)- (7.17). See [73] (appendix B), for an outline. In order 
to get (7.16)- (7.17), one performs the integrations over all the classical momenta 

) and the one over 7v c , which is also Gaussian. Then, the 



like in 



4.2 



effective classical partition function for the three-dimensional double-stranded open 
macromolecular chain, with constrained bond lengths and angles reads |73|: 



7 {ds) 
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(7.17) 



The \A&) 



N- 



-1/2 



's are given, for each s = 1, 2, by the right-hand-sides of Eqs. (4.25 ) 



(4.26) and (4.27). The function , arising from strong Va ' (hard bond angles), is 



(r) 



\d (r) u {r) + d {r) u (r) I 



sharply peaked at 

(r) (r) (r 0) (r) 

2d - d-!^- . As an approximation, we take 



d {r) 



where (dj^ +1 ) 2 



M \2 

'j+i) 



5(u 



j+ 

(r) (r) 



-1 



/3 



5 denoting Dirac's delta function. Y2 r =i ^i 

r(r) ^ 



(d)") 2 + ■ 
to be approximately proportional to 

Vd s in Eq 



(r) 



(7.17) 



is the restriction of the previous Yl r =i V\ + Vds , when all bond lengths and angles are 



constrained. Through similar arguments, Eq. (7.14) would lead to classical models for 
a ds chain, subject to stretching forces. 

Eqs. (7.16)- (7.17), with further physical approximations, have provided a basis for 



certain models of dsDNA 73, 74j: see section [8] for an outline. 



Chapter 8 



Classical effective models for 
double-stranded open DNA 

Let us consider a typical three-dimensional ds open DNA macromolecule (dsDNA) (say, 
B-DNA [2j [3] ) , at thermal equilibrium at temperature T in an interval from about room 
temperature up to about the melting or (thermal ) denaturation one, T m ~ 360 K. 
That dsDNA is formed by two open single strands of DNS (ssDNA). For T < T m , both 
ssDNA are bound to each other, forming dsDNA. For T > T m , dsDNA becomes two 
separate ssDNAs, each of which retaining still its separate existence as an extended and 
connected structure, provided that T be not too high. 

Each single strand (ssDNA) of real dsDNA is formed by a very large number N of 
nucleotides. N can vary much from one species of dsDNA to another. To fix the ideas, 
let us take: N ~ 10 10 . In turn, each nucleotide is formed by a sugar, a phosphate 
and a base (either A or C or G or T). The masses of the four bases A, C, G and T 
differ from their average mass by less than about 5, 13, 18 and 11 per cent, respectively. 
To simplify the picture, we shall regard each open ssDNA as a single discretized chain 
formed by N(~^> 1) basic units, also referred to here as nucleotides, all with equal mass 
M (such that NM equals the total mass of the DNA single strand). Thus, M includes, 
in an average or effective sense, the contributions of the masses of sugars, phosphates 
and bases. So, M is larger than the average mass of the four bases in the strict sense 
(say, of A, C, G and T). 

As stressed previously (see, for instance, [l0]), dsDNA behaves, as far as a variety 
of phenomena is concerned, as a macromolecule ressembling, as a first or rough 
approximation, a Gaussian chain in which certain effective interactions should also be 
taken into account. It is certainly most fortunate that those simplifying (Gaussian- 
like) features hold approximately for dsDNA macromolecules, which play an absolutely 
priviledged and unique role in Molecular Biology. Of course, there are other complicated 
phenomena for dsDNA, which cannot be described, even as a rough approximation, on 
the basis of an effective Gaussian chain with interactions. The latter simplifications are 
not valid for other macromolecules. 
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monomers 




Figure 8.1: Monomers in the rth-strand of a double-stranded open 
macromolecular chain. 



Based upon Eqs. 



(7.16) and (7.17) (with the i* 1 - s corresponding to hard bond 



angles) for the classical partition function of three-dimensional open dsDNA at thermal 
equilibrium, we have performed a detailed analysis and further approximations, in 
T < T m and as T approaches T m from below 73, 74 . Below, we shall outline the 



main outcomes of those studies, which seem to be supported, up to certain extent, by 
previous wisdom in [To] (and try to implement it). 

Let T < T m . Covalent forces ( the conterpart of U(y) in [3]) constrain approximately 
constant bond lengths and bond angles. In each single ssDNA in dsDNA, all bond 
lengths are taken as approximately equal to one another ( d), as are the cosines of all 
bond angles ( denoted as P^). Specifically for B-DNA, the interactions Ylr=l ^ anc ^ 



Sr=l U ^ ( section 



7) yield d ~ 0.7 nanometers and f3^ ~ 0.8, respectively. 



In each ssDNA of dsDNA at thermal equilibrium, also for T < T m , we have carried 
out medium and large distance Gaussian-like approximations, that lead to certain 
effective monomers (e-monomers), as natural molecular blocks for medium and large 
length scales 73, |74j. The analysis and Gaussian- like approximations based upon 
Eqs. (7.16) and (7.17) make those e-monomers to appear explicitly. That enables to 
regard each ssDNA in dsDNA as formed by L(= (N — l)/n e ^ N/n e 3> 1) effective (e- 
)monomers. Each e-monomer is formed by n e (> 1) nucleotides and has effective length 
d e . Lower limits are n e = 20 and d e ~ 10 nm, while upper limits of d e and n e could 



be relatively close to the persistence length dd s (— 50 nm) of dsDNA 73 74 . Different 
e-monomers behave as statistically independent from one another, in some approximate 
way ( exce pt for some weak effective residual interaction V among them, to be given in 



Eq. 



1) and discussed below). In the r-th chain, r = 1,2. 



(r) (r) 

and Rj +1 will denote 

the three-dimensional position vectors of the origin and the end of the l-th e-monomer, 

(r) 

I = 1, ...,L: no confusion should arise between the actual R, s for monomers and the 
R 4 W 's (i 



1, ...N) employed in 7.1 for the position vectors of atoms. 



All e-monomers in dsDNA are subject to effective or residual intra-chain and inter- 
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chain interactions (weaker than all covalent ones), all of which are described by an 
effective potential V ( Eq. 



distance approximations 



73 



(8.1) ), in the domain of validity of Gaussian and long- 



er) 

V depends on all . 



On the other hand, in real dsDNA below T m , residual ( intra-chain and inter-chain) 
interactions do operate and lead to longer (super-)monomers. The sizes of the super- 
monomers in the double-stranded system would be characterized by the persistent length 
dds (about 50 nanometers, amounting to 150 nucleotides, for B-DNA) or, equivalently, 
by the Kuhn length (about 100 nanometers, for B-DNA) [3|[5|[l0]. The scales of d e and 
dds are similar, although d e is certainly smaller than dd s - A ds- monomer has a length 
equal to the persistence length dds and is formed by two parallel substrands of ssDNAs. 
We emphasize that a ds-monomer should not to be confused with the single e-monomers 
in each ssDNA (with length d e ). The fact that dds > d e can be understood as arising 
from some effective interactions included in V, which would give rise to some repulsion 
in the double-stranded structure. In this connection, we remind that the interactions 
between phosphates in DNA are repulsive. For further comments, see subsection 3.4 in 
(73). 

Typical potential energies of the A — T pair differ appreciably from those for the 



C — G pair 93 . Specifically, the potential energiy of a A-T pair is smaller, in absolute 
value, than that of a C-G pair. However, the influence of those inhomogeneities in the 
effective residual interaction V among the effective monomers in dsDNA could possibly 
be rather weak and be smoothed out, due to the averaging over the interactions of the 
nucleotides included in each monomer, and we shall disregard such inhomogeneities. 
Thus, in our main analysis, dsDNA is regarded as homogeneous. 

The total effective residual intra-chain and inter-chain potential V among all e- 
monomers in dsDNA is taken as: 



V = Vo + J2v$ + V 2 . (8.1) 



Vo is the potential between all pairs of (complementary or mate) e-monomers at the 
same positions in the different strands (different ssDNAs). ^r=i e * s ^ ne potential 
between different e-monomers in the same ss DNA . It could be regarded as the result of 
averaging Ylr=i over monomers (recall ( |t. 1 )). It also takes into account excluded- 
volume effects. V2 is the potential between pairs of e-monomers at unequal or non- 
complementary positions in the different ssDNAs. 

Vo, J2r=i V\e an d ^2 are repulsive for short distances. Yll=i ^1 e an< ^ ^2 take care 
of stacking interactions. Vo and, eventually, V2 are attractive at intermediate and large 
distances. 



As a result of a series of approximations, Eqs. (7.16) and (7.17) yield the 



following effective classical partition function for three-dimensional open dsDNA at 
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thermal equilibrium, in terms of configurations of e-monomers, for T < T m 73 74 : 



Z 2 
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4^jj 



G(L). 



5.2) 



$.3) 



Mtot = 2MN is the total mass of dsDNA. The region in which the ds open chain moves 
is a sphere of very large radius Ro (compare with (4.35)). An additional factor (47ri?Q)/3 
has now been included, in order to facilitate the comparison with 74 . Such a factor, 



arising from the overall CM degrees of freedom of the whole dsDNA, has been factored 



out in all previous calculations, up to and including those in 7.1 Z^ app (one per ssDNA) 



has appeared previously in Eq. (4.34). Z# i0f , p 's are T-independent, and they will either 
cancel out or not be relevant here. Anyway, Zr app is given in Eqs. (C.4) and (C.2) in 
[73]. For the actual open ds system: 



G(L) - G(R£ fl ,R^ +1 ; R^\ R^ 2) ; L) 



.r=l V'=2 



3dW 

V 



W, 



eq 



2 L 



.r=l 1=1 



exp[-(A: B r)- 1 y] . 



W eq (8.4) 
(8.5) 



On the other hand, Wg(R^ 



7+1 



RS r) ;2^) 



denotes the Gaussian distribution for the Z-th 



monomer m the r-th strand, namely, [3/(2ttc^)] 3 / 2 exp[-3(R, ( J 1 - R^) 2 /(2dl)). Notice 



that Z and G(L) are manifestly rotationally invariant. Eqs. (8.4) and (8.5) generalize 



not only (4.35) and (4.36) but also other one-dimensional models of dsDNA 89, 91 



Through similar arguments, Eqs. ( 7.14 ) , ( 8.4 ) and ( 8.5 ) would lead to classical effective 



models for a ds chain, with e-monomers at the ends subject to stretching forces. That 
could possibly provide some basis for comparisons with experimental results observed 



in DNA stretching experiments [94, 95 and their interpretation by means of the worm- 
like chain model 96, 97 . There appears to be overall consistency between values of 
persistence length obtained in [97] with those in [5j [lOj, as it should be. A closer 
analysis lies outside our scope here. 

We remark that, for T < T m , Z and, hence, Z2 include, in principle, not only 
contributions from the bound ds structure (denoted with the superscript bo) but also 
from unbound configurations, with the two separate single strands unbound from 
each other (denoted with the superscript ub). Accordingly, it is reasonable to write 
Z = Z(H + Z^ and, hence, Z 2 = + zf\ where Z^ (zf o) ) and Z^ (Z^) 
are the contributions due to bo and ub configurations. Physically, Z^ should dominate 
over Z^ uh > the more, the larger T m — T(> 0). As T m — T — > 0, one may expect that 
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such a dominance of Z^ bo ^ disappears. The decomposition Z2 = + Z^* and those 
features are not transparent, a priori, in (8.2), but they can be supported by other 
calculations. 

Above thermal denaturation (T > T m ), dsDNA, having melt, becomes two separate 
ssDNAs. Then, for T > T m , the effective residual interactions contained in V among 
monomers belonging to different ssDNAs become negligible (as both ssDNA's are far 
from each other). Then, for both Vq and V2 in Eq. (8.1 ), one has: (KbT)^ 1 (V0+V2) — 0. 
For further discussion of ssDNA for T > T m , see [74] . We have also studied the dynamics 
of open dsDNA close to thermal denaturation, based upon both the approximations 
summarized above in this section and a Smoluchowski equation. Such a study lies 
outside the scope of the present paper, devoted to constrained macromolecular chains 
at equilibrium and, hence, will be omitted. See |74j. 
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Chapter 9 

Concluding comments 



We have concentrated on quantizing Fraenkel's model (very stiff flexible chain), by taking 
54, 6l] and, mostly, [56] : see the last paragraph of subsection |2.2l In a 



into account 

nutshell: upon quantization, does one obtain the physically expected large quantum 
zero-point energies of hard degrees of freedom? Do those zero-point energies depend 
on soft degrees of freedom? To the best of our knowledge, no exact or asymptotically 
exact answer, operating directly on the Schrodinger equation and providing affirmative 
answers to those questions, exists at present. However, a moderate step forward has 
been made in the last years |26l 67 71-75, 86 for equilibrium quantum-mechanical 



partition functions of three-dimensional macromolecules, based always on the same 
underlying physical ideas. Specifically, upon quantizing Fraenkel's model and applying 
Peierls' quantum- mechanical variational inequality for large vibrational frequencies, the 
expected and physically correct large quantum zero-point energies of the constrained 
(hard) degrees of freedom are obtained exactly. They are proved analytically to be 
constant (and, hence, independent on the unconstrained soft variables) for several 
ss macromolecules (open, closed-ring, freely-jointed, freely-rotating, etc) and for ds 
open ones. We stress the crucial importance of choosing the variational trial wave 
functions adequately: otherwise, one may obtain either contributions depending on 
the unconstrained soft variables or results which, even if independent on the latter 
variables, do not coincide with the physically expected quantum zero-point energies (see 
the comment after (6.11 )). The proof of constancy becomes considerably complicated, as 
one proceeds from the open freely-jointed chain to chains with constrained bond lengths 
and bond angles or closed-ring constraints: we have provided in appendices A and D 
some (so far, unpublished) details which complement, in outline, various aspects given 
previously (71} 73 , [75] . That constancy leads to specific quantum partition functions 
and hamiltonians (Zq and Hq) for the soft angular variables, neatly separated from the 
hard degrees of freedom. One always finds structures similar to those in Eqs. (6.9), 
(6.10) and (6.11), with different Wang (associated to all constraints except those for the 
constrained bond lengths): see section [6] 

We are aware of the increasing complexities of the Zq , s and Hq's derived for the 
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successively quantized chains, with increasingly complicated constraints, in sections [4] 
through [7] Our justification for them is twofold. First, Peierls' inequality yields, after 
delicate exact cancellations, the correct constant quantum zero-point energies. Second, 
the resulting Zq's and Hq's are rotationally invariant and the total quantum angular 
momentum is conserved, which are very important physical properties. In view of 
those nontrivial results, the quantum variation inequality appears to operate in the 
right direction even if, recognizedly, it may not be providing us, as yet, the definitive 
quantum- mechanical formulation for constrained macromolecules. 

What practical consequences do Zq's and Hq 7 s imply? How do their applications 
differ from and/or improve those from cCHDa? For that purpose, by taking the classical 
limit, the Zq's become the classical partition functions Zq's, which are, respectively, 
different from the classical partition functions found in cCHDa for similar chains. The 
Zc's for open chains yield, after certain large-distance approximations, several quantities 
(bond-bond correlations, squared end-to-end distance, probability distribution for the 
end-to-end vector, and others), which agree consistently with those from the standard 
Gaussian model in Polymer Science. Since a similar consistency is met with cCHDa, 
the comparison with the standard Gaussian model does not distinguish, thus far, the 
consequences of our variational QMa from those from cCHDa. 

As another methodological application, we have also applied the above 
Peierls'inequality approach plus classical and long-distance approximations to three- 
dimensional ds macromolecules, like dsDNA. We have derived the corresponding Zc's, 



bearing structures which generalize other one-dimensional models of dsDNA [89 91 
One could also develop similar quantum- mechanical approaches, starting from Eq. (7.1) 
with similar Y2r=i > Ylr=i Ua (so as to constrain bond lengths and angles) but 
with some suitably chosen Ylr=i ^i +^ds ; so as to generate other kinds of (eventually 
weaker) angular constraints in each strand among neighbouring atoms (which be neither 
nearest-neighbours nor next-to- nearest-neighbours) . Then, by extending the approaches 
in sections [7] and [8j one could get eventually other models for dsDNA. 

Thus, whether our variational QMa has practical applications for macromolecules at 
thermal equilibrium which differ neatly from and/or improve those from cCHDa remains 
open, thus far. 

We have limited ourselves to static properties (equilibrium partition functions) upon 
quantizing very stiff flexible chains. We have not undertaken the far more difficult task 
of quantizing very stiff flexible macromolecules off-equilibrium (that is, of analyzing 
dynamical properties ab initio, at the quantum level). Anyway, a model for the non- 
equilibrium evolution of a ds macromolecule, based upon the Smoluchowski equation 
and the approximate Z obtained in [73], has been proposed in [74]. In so doing, we have 
followed a pragmatic procedure (not to be confused with a first-principles treatment): 
having obtained an approximate classical partition function (Z), we have constructed 
directly the classical Smoluchowski equation which has as equilibrium solution the 
distribution function characterizing uniquely Z. This model has been employed to 
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study thermal denaturation of dsDNA: specifically, an approximate formula for the time 
duration required for thermal denaturation to occur (about the melting temperature) 
has been obtained. That pragmatic method appears to be consistent, at least in spirit, 
with that followed in other approaches to the dynamics of constrained macromolecules: 
compare, for instance, with [7, 18 . 



Biological macromolecules and the processes which occur in or involve them give rise 
to a fantastic variety of phenomena. A good number of the latter may well be described, 
to a sufficient degree of approximation, by classical approaches (cCHDa). At certain 
stages in the analysis of macromolecules, it may be unclear whether quantum-mechanical 
approaches have enough practical consequences which differ from and/or improve those 
from classical ones. But it is impossible to accept that classical formulations will account 
for all the above phenomena, at all scales down to the nanometer one. Thus, and as 
a matter of principle and of scientific strategy, it seems that, if possible, one should 
pursue on disposing of both classical and quantum formulations. In fact, at some short 
spatial and/or temporal scales (at some adequately large energy scale), one will have 
to deal with some genuinely quantum features and to analyze them. See the comments 



in sectiorj3l in connection with Schrodinger's book 69 . This has been a leitmotif to 



motivate the researches reported in this tutorial review. 
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Appendix A 



Single-stranded open freely-jointed chain: proof of (4.13)- 

(P3 



We shall use [dy] = [n^ 1 yjdyi] [dO] and [dO], *(y) and <Myz) ( E q«- ([^13]), ((45J) 
and (4.7), respectively). We shall consider successively all contributions from: 



JV-l 



JV-l 



JV-l 



i=l 



i=l 



i=l 



Hq^u being given in Eq. (3.5). First, we shall suppose that U(y) is the sum of Morse 
potentials in Eq. (4.3). We shall start with the simplest structures. The contributions 
from 2~ 1 Bi{&i -&i/y 2 ] and a^-i ■ &i / {Mi{yi-\yi)) follow immediately, by letting ui — > +00 
and using Eqs. (4.10) and ( 4.11[ ). One gets: 



[dii\Mo><p) 



/ 



(A.2) 
(A.3) 



Mi di-xdi 

Next, the contributions from M~ [h 2 Ui-i-Ui(d 2 /dyi-idyi) — ih((ui-i-ai)/yi)(d/dyi-x) — 
i/&((a*_i • Ui)/yi-i)(d/dyi)] follow by integrating by parts over y^ and yi-i, letting 
ui — > +00 and using Eqs. (4.10) and (4.11). One finds: 



1 f /i 2 Uj_i • Uj iHui—i ■ &i ihsLi—i ■ Ui 



Mi 



^{6, tp) . (A.4) 



di-\di di_\di di-\di 

We shall now deal with somewhat more complicated structures. We shall make use of 
the following identity for an arbitrary radial function <p(y): 

d<j>(y) = 1 dct>{y) 2 
2 8y 

We consider: 



(A.5) 



JV-l 



JV-l 



Hyfdy^dmUMyi^Ao^)} 



1=1 



1=1 



2 



1 8 N ~ 1 
•2SV^(nM!/l)i(^)l(A.6) 



y% dyi 



i=i 
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We apply Eq. (A.5), integrate by parts over yi, notice that the contributions from y 
and from yi 
result is: 







+00 vanish, let w/ — > +00 and use Eqs. (4.10), (4.11) and (4.12). The 

Bih 2 



2d 2 



(A.7) 



We treat, also as wj 
N-i 



+00, the remaining contributions: 



JV-1 



AT-1 



JV-l 



I] y fd yi ][dn][ll Myi)^]*lYl { -^~^ + u(y)\ U UvM (A. 



i=l 



1=1 



i=l 



There are more than one way to treat them. See |26|. Possibly, the simplest one is: 



- (-Bi)h 2 d 2 

2 dy 2 



+ V M (yi)}4>M,i,n=o(yi) -> E M ,i,0^>M,i,n=o(.yi 



(A.9) 



EM,i,n=o is given by (4.9 ), with the corresponding D{, Oj, Wj (see the comment justifying 



(4.9)). By using Eqs. (A.9), (4.10) and (4.11) in Eqs. (A. 8), by recalling the results in 
Eqs. (A. 2), (A. 3), (A. 4), (A. 6), (A.7) and by defining the new variable e$ = ih\ii — aj, 



we arrive at Eqs. (4.13), (4.14) and (4.15) for the single-stranded open freely-jointed 
quantum chain. In so doing, use is made of: Uj • aj = and of aj • Uj = 2ih. 

Next, we shall turn to the case in which U(y) is the sum of harmonic-oscillator-like 
potentials given in Eq. (|4.1|). As uji — > +00, one has: 



A-Bi)h 2 d 



i-07 



dy 2 + 2 Bl {m 



di) 2 ](pi(yi 



IjjjJ; 



(A.10) 



The results in Eqs. (A. 2), (A. 3 and (A.7) continue to hold. The main difference now is 



that Y,i=i E M> i ;n=0 is replaced by £-=i Hwi/2. All that yields (|4.13)- (4.14) 



We remark that neither ihu[ nor a; are Hermitean operators but, on the contrary, 
= ih\i[ — en is. Using Cartesian components, we write e/ = (ej i, e^2, e^s) and 
u l = (fj 1, Wj2,uj 3). Then, one can derive the following "angular commutation 
relations": [ui )CC , ejk j( g] = ih5ik(Sa/3 — u^aUifl) 86 . It suggests that Cartesian components 



of u; could be regarded as conjugate variables (quantum-machanically) of e;, in some 
extended sense. Those (purely algebraic) "angular commutation relations" ( together 



with the closed algebra formed by (4.16) and the standard commutation relations for 
orbital angular momentum in subsection 4.1) hold independently of the uji — > +00 



limit, and do not look particularly interesting before taking the latter limit, because 
radial variables also matter. However, after w; — > +00, all those algebraic relations 
constitute a distinguishing feature of the remaining ( unconstrained) angular variables. 
Whether or up to what extent they constitute strict quantum equivalents of classical 
radial constraints constitute open questions. 



Appendix B 



Weak next-to-nearest-neighbours interaction: G(q) and 
<(x» - Xl ) 2 ) 



For small q, the dominant contribution to G(q), as given by the ratio in Eq. (4.38), can 
be evaluated approximately as follows. Studies of [A^r-i] -1 / 2 for large N have been 
carried out, with successive improvements, in [26] (Subsection 5.4), [72] (Appendix B) 
and [73] (Appendix C): in short, [Atv-i] -1 ^ 2 takes on its dominant contributions when 
all (uj 



u, 



,iV-l 



= 1, i,j = 1, ..N — j. Then, in each integral over 9 S , s = 2, 

in both the numerator and denominator in Eq. (4.38), we keep only the contributions 
over two adequately small intervals of size 59 s about 9 s -i and it — Q s —\. So, we extend 



the arguments in 26, 72 73 , with two modifications: i) we impose the classical limit 
restrictions, by including a factor pq for 9 s+ \ ~ ir — 9 s , ii) the contributions of the values 
taken by U nnn when 9 S is close to either 9 s -\ or tt — 9 s -\. The only contributing values 



of v are v(2d) (for u» • Uj+i ~ +1), and v(Xth/ '{2 l / 2 d)) (for u» 



Ui+i ~ -1 + \ th 



/2d 2 



Thus, having integrated over all 9 S , s = 2, . . . , N — 1 (by retaining only two adequately 
small regions in those integrations), one gets the following approximate representation: 



G(q) ~ 



Jq sin9 1 d9 1 z(9 1 ;q) N 
Jq s'm9id9iz(9i;q = 0) N 



(B.l) 



z (Qi \ q)n turns out to be a sum of 2 N 2 terms: 



)N 



exp [(N - 2)(k B T)~ 1 v(2d)] 



N-l 

n 

.8=2 



^2a i l N \p)exp[iqdlcos9 1 ] . (B.2) 



J l 



p was given in subsection 
1)(P) : 



below [with 



4.3 



a i (p)(> 0) are certain polynomials in p characterized 
0]. The summation in Eq. (B.2) should be understood 
as follows: for N even (odd), I = 2k + 1 (I = 2k), the integer k ranging from 
-(N/2) + 1 up to (N/2) - 1 (from —(N — l)/2 up to (N - l)/2). The aj^O^'s, 



following from the above integrations, are obtained successively from the following 



recurrence: a 



(N+l), 



(p) = a^l(p) + pa^J/_ 1} (p) 



(TV) 



(N+l) 
(7V-2V 



(p) = pa^\(p). For instance, a^'(p) = 1 



with a 

(4) 



(N+l) 



N 
(P) 



(P) 
P + P 2 



a 



SJ_iG°)> and 

a { f[(p) = p. 
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Eqs. (TO) and (Bib yield 



(JV ), ,sin[(iV-l-2j>d] 



(AT - 1 - 2j)$d 



/TV) 



(AT) 

N-1-2J 



(p) 



(JV) 

-(AT-l-2i) 



(p) 



(1 + P) 



iV-2 



(B.3) 
(B.4) 



The summation in Eq. (B.3) should be understood as follows: for N even (odd), the 



integer j ranges from up to (N—2)/2 (from up to (N— 1)/2). The following properties 



hold: EjPf\p) = l,/3f J (l) 



if N is odd, except /3(^ 1 )/ 2 (1 



(N-l 



(7V-l)/2^ 



)/2 

JV-l 
i 



7V-2 



if TV is even, and ^ N) (l) = ( 7 3 
)/(2 N ~ 1 ). Notice that, upon performing the 



fN-l 



)/2 



JV-2 



inverse Fourier transform of Eq. (B.3), W(r) ~ (47rd | r |) 1 Y^jj (p)^(l r I ~jd) 



(5 denoting Dirac's delta function). This representation for W(r) strongly reminds 
(although it does not coincide with) the end-to-end distribution appearing in the 
random-flight approach to polymers (6j. Eq. (B.3) yields: 



((xjv - xi) 



Using the properties of the (3j N \p), one finds the following expressions for the case N 



d 2 T] 



N-2 



(1 + P) 



N-2 ' 



c) V 



(B.5) 



3=0 



even (the case N = odd being analogous): = [(N — 1 — j)(N — 2 — j)(iV - 
2)j!(iY-l-j)!] _1 if j is odd, while cf ] = [(N -1- jf{N -l)!][(i + l)j\(N - 



1)!][(J + 



1-J 



ill- 1 



if j is even. The computation of rj, leading from ( |B.5 ) to (4.39), is outlined in appendix 
C. 

We claim that the approximations leading from Eq. (4.38) to Eq. (B.3) hold 



for small q, say, up to and including order q (and, so, they are valid to compute 
((xjv — xi) 2 )). This statement is supported by the consistency obtained with the 
standard Gaussian model up to order q 2 when U nnn = [6] . 
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Additional computations for subsection |4.3| and Appendix 
B 



Without loss of generality, we limit ourselves to the case N even. Then, some direct 
algebra yields 



r, = (jv-i) 2 [i + ri] + (jv-i)jvr 2 + r 3 - (jv-i)t 4 

- 2- x [2(iV- 1) + 1]T 5 . 



(C.l) 



The quantities Tj, i = 1,...,5 are certain sums involving ratios of factorials 
(combinatorial coefficients). We shall give below their initial and final expressions. 



1 



(N-iy. 



2t + l(JV-l-2t)!(2t)!' 



{N-2)/2 

E 

1=1 

(N-2)/2 

2i + l(JV-2i)!(2t-l)! 

(l + p )AT + (1 _ p) iV p N-l 



„2i 



(1 + pf-j^- p) 

2Np 



N 



I 



1 



(1 + P) 



AT+l 



(1-p) 



AT+l 



2iV> JV + 1 

Ti + T 2 + (TV - l)p [(1 + p) 
(1 + p)^- 1 + (1 - p) N ~ l 



2N(N + l)p 2 

_ 2i (N-l)\ 2i _ x 

(N - I - 2i)\{2i)\ P + {N - 2i)\{2i - l)\ P 

JV-2 p N ~ 2 ~\ 



+ 1 , 



(N- l)\p 



l/i 2 *- 1 



(N-2)/2 

E 

i=l 

[(l + p) 7V - 1 + (l-p) Ar - 1 ]-2-2T 1 , 

(JV-2)/2 



2i + l(JV-2»)!(2i-l)! 



4i 



5^ 2/ + I ! A' - i 



(JV- l)!p 2i 



[(1 + P) 



AT-l 



(1-P) 



AT-li 



2p N ~ 1 - 2T 2 . 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



(C.6) 
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When performing those sums, use has been made of the following formulas 

j.l , h 4- n\ m + l - 1 



and 



y— 



>>>■ _ p k = i l + P) 



k=0 



(m — k)\k\ 



\ - ml k 

f^k + 1 (m-k)\kV 



(m + l)p 

1 - (1 - p) m+1 
(m + l)p 



(jV-2)/2 



E 



(^-1)1 p 2i-l 



' 2»(2» + l)(iV r -2*)!(2f-l)! 



1 



(iV + l)iVp 2 

(AT-2)/2 



7V+1 



(1-P) 



E 



2/ 



(7V-1)! 



(N - 1 - 2t)!(2i)! 



p* + (2i - 1) 



(JV + lJp-^ 1 
(N-l)\ «_! 



i=i 

:(iV-l) P [(l+p) 



(JV-2»)!(2t - 1)! 



(C.7) 
(C.8) 



(C.9) 



V-2 _ p^- 2 ! 



(CIO) 



Upon combining Eqs. ( C.2[ ) through Eq. (C.6) and further substitution in Eq. (C.l ), 
one finally arrives at Eq. (4.39). 
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Single-stranded closed-ring freely-jointed chain: proof of 



(5.5) 



In this appendix, we report some details and the results of the evaluation of (Hq) 



(HQ,in) + (U(y))- Use will be made of harmonic oscillators (Eq. (4.1)), of the unit vector 
uat defined through: = — X^i" 1 Yi = Un^-n and of the methods in Appendix A. The 



results given below will come from Eqs. (3.5) and (4.1). Using Successive integrations 



by parts yield (as all frequencies Ui, i = 1, . . . , N, grow very large): 



dy 2 



1=1 



i=l 



2dj 



N-1 

+ E 

i=i 



~Jb n ~ 



'dft]\^(e,ip)\ 2 (ui-u N 



2 S(yN - d N ) \ 

! 4, I 



E 

i=i 

N-1 

E 

i=i 



2d t 
h 2 Bi 



d_ 
dyi 

dyf 



'dfl]\M0,<P)\ 



,S(y N (0,cp;yi) - d N )\ 



d 2 
a N 



IdCl] 



7 



Vi=di 



(D.l) 



E s-» 2l #- 



E 

1=1 



2d 2 ' 



(D.2) 



N-1 

E« 

i=l 



jV-1 

+ V — 

«=i * 

iV-l 

«=i 1 



;dn]V;(e,^) 8^,(0, ¥>)] 



5{vn - d N ) 
d 2 

a N 

Sjm - d N ) 
d 2 

a N 
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JV-l 



+ E 



i=i 

N-l 



+ E 



i=i 

N-l 



4£?jv 
±B N 



J[dn]\M0,<p)\ 2 (vo i -"N)'- 



d 2 

a N 



2 5(y N - d N ) | 



d 2 
a N 



E 

i=2 
JV-l 

-E 

i=2 
JV-l 

+ E 

i=2 

TV— 1 



ft 2 



Mi 

Hun 

2M, t B N 

h 2 

4Mi 



Uj-l • Ui 



i=l * 

a 2 



Mvn ~ d N ) 
d 2 

a N 



, (D.3) 



dyi-idyi 



J [dft]\i/j a (9, ^)| 2 (u;_i • Ui)(uj_i • UAr)(ui • u N ) 6<yVN d2 ^ j 



d 2 



dyi-idyi 



[dnn^^^l^ui-i-ui) 



S(y N (e,(p;y i - 1 ,y i ) - d N ) 



d 2 

a N 



Vi=di 
Vi-i=di-i 

(D.4) 
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-E 
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N-l 



-E 
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JV-l 



1 a^i • 

Mi yi-iyi 



2M t B N 



i=2 
JV-l 



£ 2^ {/ • u iV )(u, i _ 1 • u e JK • ujv) ^^ } 



=2 

JV-l 

E 

i=2 
JV-l 



HUN 

2M { B N 



j [<m]\i>a(6,<P)\ 2 (u Vi -i •U/v)(u^_ 1 ■U Vi )(u Vj 

E to {/ i dS! i«< s - ^ [<*-■ ' 



d 2 
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x%JV - 




} d 2 

a N 


djv) J 


J(y N - 




' d 2 

a N 




J(y N - 





i=2 
JV-l 



-E — 



i=2 
JV-l 



[&n\\M9,<p)? 



(aj-i • aj 



<%jv - djv) 
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-Esm^I/^*''''')!'''-''''-^^- 



i=2 
JV-l 



<%jv - <jjy) 
; d 2 

a N 



ai_i 



<%/v - djv) 
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(D.5) 
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N-l 
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N-l 



Mi 
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1 Vi dy 



aj_i • u, 



i-i 



j a_ 

2/i-i % 



2MB 



i-DiV 



[clO] | -00- (6> , | 2 (Xli—JL • Ujv)(Uj_l • U e .)(u 0i • Utv 



S(vn - d N ) 



d 2 

N 



2M-B N { / l dn ]\^( 9 ' ( fi)\ 2 ( u i^ • u N){ui-i ■ u w )(u^ • ujv) ^^^ ^ j 



i=2 
AT-1 



-E 

i=2 
NX 

-E 

i=2 
AT-1 

-E 

i=2 
N-l 

-E 

i=2 
AT-1 

+ E 

i=2 
AT-1 

+ E 

i=2 
AT-1 

+ E 

i=2 



2MiB N 

Hun 

2M~B~n 
ih 



'dn}\ip a (9,ip)\ 2 {Ui ■ u N )(ui ■ U0 i _ 1 )(u e ^ 1 ■ u N 



^ 5{y N - d N ) \ 



d 2 

N 



dil]\ip a (0, (p)\ 2 (ui ■ u N )(ui ■ Uy^^Utp.^ ■ u N ) 



S(yN ~ d N ) \ 
d% 1 



4Mi 

ih 
4M. 



d 

dyi-i 
d_ 



[dnuAQ, ^)l 2 { —, ' a ° 5{m{6, ^ y r l) " dN) 

di dpq 



Di-l=di 



dfi]|Vv(^)| 



2 (uj • aj_i) 5(y N (6, 99; jfe) - d N ) \ 



di-i 



d 2 

a N 
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Mididi_i 
ih 

Mididi-i 

ih 

2Mj 



dn]|^(e,^)| : 



(Ui-i ■ a.,) (u, • aj-i) 
didi-i 
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, (D.6) 



AT 

E 



2Bi {Vi 



di 



N 

E 



(D.7) 



Besides g/jv (the meaning of which is given in Eq. (5.4)), the above equations contain two 
new functions, namely, Un (&■><£'■> Hi) and Vn(G, <P', Vi-iiVi), which depend on the variables 
displayed respectively. In turn, those two functions are defined as follows: 

for all yh = dh, h = 1, .., i — l,i + 1, .., N — 1, 



that is, h 



v^AT-1 

luLj=x yi ■ yj 



i is excluded, 



1/2 



b) y N (8,ip;yi-i,yi 



v^AT-l 



1/2 



!y/ •yjl , for all = d h ,h = 1, .., i-2, £+1, .., JV-l, 
that is, both /i = j — 1 and /i = £ are excluded. 

The meaning of, say, [a.iip a (9, (p)] [ai5(yiy — djy)] in Eq. (D.3) is the following. First, 
the differential operator aj acts upon ip a by regarding all angular variables 9i, <pi, 
I = 1, .., N — 1, as if they were independent on one another (that is, as if the constraint 
fi(UN — d^ were not operative). On the other hand, the interpretation of a.i5(yN — djv) 
is similar to that of the first derivative of S(f(x)) with respect to x (f(x) being a given 
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function of x): d5(f(x))/dx = (df(x)/dx)dS(f)/df, d8(f)/df being the first derivative 
of 5(f) [65j [66]. Then, aj<5(yjv ~~ <^v) a l so embodies the essentials of the closed-ring 
constraint. After the operator a^ has acted upon ip a with the understanding explained 
above, then a.i5(yM — d^) acts and implies that there is one relationship among those 
2(N — 1) angles. A similar interpretation applies for other therms containing differential 
operators with respect to angles acting upon 5(y^ — djy) and upon ip a . 

Eqs. (D.l) and (D.7) display the simplifying feature that all frequencies cjj, 
i = 1, . . . ,N— 1 go multiplied by constant (angle-independent) factors. However, in Eqs. 



(D.l), (D.3)- (D.6) one sees terms linear in con, which go multiplied by complicated, 
in principle angle-dependent, expressions. We shall prove that there are two crucial 
cancellations, when we add all those matrix elements which contain the frequency 
u>at times various functions (displaying angular dependences). The first cancellation 
is obtained if we add the corresponding terms of Eq. (D.l) and Eq. (D.3), which 
contain ojn times the following angle-dependent factors: 
N-l 



E 



dfi]|W(^)| 2 [(Ui • U7V ) 2 + (ity • u N ) 2 + (u w • u N f 



5(y N - d N ) 



d 2 

a N 



N-1 

E 



fkjJN Bi 

4B N 



d 2 
a N 



N-1 

E 

i=l 



huj]\} Bi 
±B N 



(D. 



as u. 



ifr — 1. The second (and rather non-trivial!) cancellation arises if we add the 
corresponding terms of Eq. (D.4), Eq. ( |D.5 ) and Eq. (D.6), which contain ujn times 
the following functions containing angular dependences: 
N-l 

[dft]\'ip a (9,(p)\ 2 [(ui_i • Uj)(ui_i • ujv)(uj • u N ) 



hujpf 



i=2 1 JV 

+(u 0l _ 1 • ujv)(u e . 
+(u Vi _ 1 • Ui\r)(iv 

+ (Uj_i • UAr)(ui_ 



U 0i )(u ei • U N ) + (ug i _ l ■ U iV )(u 0i _ 1 

■ u ei)(u ei ■ U N ) + (u^_ x • ujvXu^ 



u 



UN) 
■ U N ) 



■ VL0.)(VL0. ■ U N ) + (Uj_i • Ujv)(ui_i • U^Jfu^ • Uat) 
+ (Uj • U N )(\li • U 6 » i _ 1 )(U0 j _ 1 • U N ) + (llj • UAr)(uj • U w _ 1 ) (tl^.^ ■ Ujv)] 



SjjjN ~ d N ) j 



d 2 

a N 



N-l 

E 



fojJN 

2MJ3 



N 



>5(y N - d N ) | 



d 2 

N 



N-l 

E 



2MB 



N 



(D.9) 



Both cancellations in Eqs. (D.8) and ( |D.9 ) are due to the following properties: i) the 
vectors Ui,u tpi ,ug i are orthonormalized (as are Ui_i,u ¥ , i _ 1 ,u^ i _ 1 ), ii) u 2 ^ = 1. 



Notice that if we add the first term of the right hand side of Eq. (D.l) plus those 



in the right hand sides of Eq. (D.7), Eq. (D.8) and Eq. (D.9), we find: 



iV - 1 * N * 

1=1 4 = 1 



AT-1 , AT-1 
1=1 1=2 



ftkJN 

2M,B N 



N 

E 



(D.IO) 
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which is the total vibrational zero point energy of the closed chain in its ground 
state. Use has been made of the expressions for Bi, i = 1...N. Thus, all the above 
cancellations have enabled to transform the initial very lengthy expressions (containing 
angle-dependent expressions, multiplying the frequencies) into somewhat shorter ones, in 
which all coefficients multiplying the frequencies are constant. The resulting expressions 
are still somewhat lengthy. The final result is collected in appendix A in |71| . Eq. ( 5.5 ) is 
nothing but a simplified way of presenting th e fin al expression, which displays explicitly 

only the most relevant terms: / oing(^) / in (5.5) denotes the remainder, which follows 



directly from appendix A in |71] . 

If Morse potentials are employed, instead of harmonic-oscillator-like ones, the 
computations and cancellations in this appendix continue to hold, with Huii/2 replaced 
by E M ,i,n=o, i = l~.N. 
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